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ABSTRACT 

Context. Complex molecules such as ethanol and dimethyl ether have been observed in a number of hot molecular cores and hot 
corinos. Attempts to model the molecular formation process using gas phase only models have so far been unsuccessful. 
Aims. To demonstrate that grain surface processing is a viable mechanism for complex molecule formation in these environments. 
Methods. A variable environment parameter computer model has been constructed which includes both gas and surface chemistry. 
This is used to investigate a variety of cloud collapse scenarios. 

Results. Comparison between model results and observation shows that by combining grain surface processing with gas phase chem- 
istry complex molecules can be produced in observed abundances in a number of core and corino scenarios. Differences in abundances 
are due to the initial atomic and molecular composition of the core/corino and varying collapse timescales. 

Conclusions. Grain surface processing, combined with variation of physical conditions, can be regarded as a viable method for the 
formation of complex molecules in the environment found in the vicinity of a hot core/corino and produce abundances comparable to 
those observed. 
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1. Introduction 

A number of complex molecules have been discovered in the 
interstellar medium. First was methanol, CH3OH, (Ball 1 19701 1 
followed by dimethyl ether, CH3OCH3, (Snyder et al. [T974] l and 
ethanol, CH3CH2OH, (herein denoted C2H5OH) (Zuckerman et 
al. 1975 ). A variety of other isomers, isotopic variants and sim- 
ilar compounds have also been seen. Charnley et al. ( |1995t re- 
view the distribution of complex molecules while Ikeda et al. 
( I2001l & r2002l) present more recent observations. CH3OH is seen 
in cold dark clouds with a relative abundance of 10 (with re- 
spect to H2) while in hot cores this can be as large as 10"*'. The 
highest abundances are observed in regions where grain mantles 
are likely to have recently evaporated, strongly suggesting man- 
tle processing either directly or in the formation of precursors. 

A variety of computer models have been constructed to in- 
vestigate molecule formation. Predominantly these are based 
on chemical reaction networks and use a "flat" (fixed) param- 
eter space. As with the models presented here most are "single 
point", where a single representative point in a cloud is mod- 
elled. This is satisfactory as chemistry at another point with dif- 
ferent physical parameters could be modelled simply by chang- 
ing the parameters. Over time the size and complexity of the 
reaction networks has increased as more reaction rate data is 
published in the literature and greater computer power becomes 
available to process more complex networks. (See Millar 019901 1 
for a review of model development.) Millar et al. (,1991bi) model 
the production of CH3OH by a fixed parameter gas-phase-only 
model and are unable to produce abundances in excess of 10"^ 
lending further weight to the involvement of grain mantle pro- 
cessing. C2H5OH is seen in fewer sources than CH3OH and 



all known sources are dense, hot core star-forming regions, ex- 
actly the places where grain mantles are likely to evaporate. 
Abundances are in the range 10^^-10"^ (with respect to H2). 
CH3OCH3 is often seen in the same sources as C2H5OH and 
has recently been detected in hot corinos (Ceccarelli et al. 2004| 
BottinelU et al. I2004I I. It has an abundance range of 10^^-10"^ 
and the CH3OCH3/C2H5OH abundance ratio differs greatly be- 
tween apparently similar sources. (Abundances with respect to 
H2 for selected sources are listed in Table[T]) 



Table 1. Observed fractional abundances of complex molecules 



Source - Hot Core 


CH3OH 


C2H5OH 


CH3OCH3 


NGC6334F 


2.0 X 10-' 


9.0 X 10-" 


4.0 X 10-** 


G327. 3-0.6 


1.0 X 10"' 


1.0 X 10-** 


3.0 X 10-** 


G31.41+0.31 


9.0 X 10-** 


2.0 X 10-** 


2.0 X 10"* 


G34.3+0.2 


9.0 X 10-** 


6.0 X 10-' 


1.0 X 10"** 


GlO.47+0.03 


2.0 X 10 ' 


1.0 X 10-* 


3.0 X 10"** 


Sgr B2 (N) 


2.0 X 10-** 


1.0 X 10-' 


7.0 X 10-1" 


DR21(0H) 


1.0 X 10-** 






W3(H20) 


4.0 X 10-** 






W51 el/e2 


3.0 X 10-' 


9.0 X 10-' 




Orion Compact Ridge 


2.0 X 10-' 






Orion Hot Core 


l.Ox 10-« 


2.0 X 10-** 




Source - Hot Corino 


IRAS 16293-2422 


1.0 X 10-' 




2.4 X 10-' 


NGC1333-IRAS4A 


< 1.0 X 10-** 




< 2.8 X 10-* 


NGC1333-1RAS4B 


7.0 X 10-' 




< 1.2 X 10-'^ 


NGC1333-1RAS2A 


3.0 X 10-' 




3.0 X 10"* 



Data from Ikeda et al. l2001l Ikeda et al. 120021 Bottinelli et al. l2007r 
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Even under very favorable conditions fixed parameter gas- 
phase-only chemical models produce maximum abundances for 
C2H5OH and CH 3OCH3 of the order o f 10'" (Herbst & Leung 
[19891 Millar et al. ll991al Charnley et al. lT992] l. between two and 
five orders of magnitude below that observed. Again this sug- 
gests a scenario in which mantles form on dust grains at some 
point in a cloud's lifetime when conditions are favorable. The 
mantles are active and complex molecules (or their precursors) 
form on them more efficiently than in the gas phase. Later the 
cloud evolves becoming warmer and denser and the contents 
of the mantles pass back into the gas phase. Several other gas- 
phase-only models (for example Charnley et al. 1992; Caselli et 
al. [1993 ) have demonstrated that the chemistry in hot cores is 
active, with particular gas phase reaction channels initiated with 
simple species thought to originate from mantles, a possible ex- 
planation for observed complex molecule abundances. 

A number of more in-depth models have included surface 
reaction chemistry. One of the first has been developed by 
Hasegawa et al. ( 119921) and Hasegawa & Herbst (I1993I ). This 
model has 274 chemical species with 2928 reactions and in- 
cludes 1% (by mass) dust grains in the gas cloud. Molecules 
can freeze out onto grains and later pass back into the gas phase. 
Once on a grain heavier species molecules are held stationary at 
lattice binding sites while lighter species, most commonly hy- 
drogen atoms, can migrate around the grain. When mobile light 
species encounter fixed heavy species the two can combine to 
produce progressively larger molecules. While large molecules 
can thaw off the grain surface there are a number of other pos- 
sible desorbtion mechanisms that can act as well. (See Section 
12.31 ) One additional benefit of this approach is that grain surface 
catalysis of hydrogen molecules can be modelled directly. 

Shalabiea & Greenberg (1 19951 1 then take the next logical 
step in interstellar cloud modelling by including variations in 
physical parameters. They use a gas phase reaction rate network 
and interchange between gas and surface similar to Hasegawa 
et al. (II992I 1 and Hasegawa & Herbst (il993j . However certain 
physical parameters (eg. density) can also change. Shalabiea & 
Greenberg (1995) model these by including differential equa- 
tions for them and solving for the relevant parameter(s) at each 
point in time in much the same way as species abundances are 
solved for. They also coin the terms "pseudo-time-dependent", 
where the model chemistry evolves over time while the param- 
eters are static and "real-time-dependent" where both chemistry 
and physical parameters evolve. Their paper provides compar- 
ison between the two types, with the gas and grain model in- 
cluding 218 chemical species and 2075 reactions. With the more 
recent identification of "hot corinos" (Table [1] & Section 13.1b 
where the temperature range is optimal for the grain surface pro- 
duction of larger molecules and their subsequent return to the 
gas phase some specific modelling of these type of objects has 
been done. Garrod & Herbst (12006!) use a gas and grain reac- 
tion network with 655 species and 6509 reactions. The same 
model was used by Aikawa et al. (I2008I I. Subsequently Garrod 
et al. (I2OO8) produced an extended and more generalized ver- 
sion. Comparison of the output of these models with the results 
presented here is made in Section|5] 

As is discussed further in Section 12.31 some mechanism(s) 
must be returning material from grain mantles back to the gas 
phase since the alternative is depletion of molecular species 
heavier than H2 in a timescale shorter than the lifetime of known 
clouds, which is not seen. More recent modelling work has fo- 
cused on the investigating the effectiveness of several proposed 
desorption mechanisms. Willacy & Millar ( 119981 ) provide a set 
of gas and surface models which include different desorption 



mechanisms and compare their effectiveness. These models use 
282 species and 4864 reactions. Garrod et al. (120071 ) specifically 
model formation of CH3OH in a quiescent cloud and point out 
that if grain surface processes are invoked to explain the abun- 
dance of CH3OH then at least one non-thermal desorbtion mech- 
anism must be active. Tielens & Charnley (119971 ) have argued 
that, at least in some cases, the reaction network model can- 
not be justified and that a Monte Carlo simulation approach is 
more appropriate. This has led to some debate in the current 
literature about the relative merits of the two methods. Willacy 
& Millar (119981 ) discuss this and conclude that currently there 
are computational impracticalities in using the Monte Carlo ap- 
proach in real-time-dependent models, and further that the two 
different approaches may give fairly similar results for systems 
of increasing complexity. A direct comparison of the two ap- 
proaches by Garrod et al. ( 120091 ) demonstrates considerable sim- 
ilarity between them. Certainly previous gas phase only reaction 
rate models have yielded results comparable with observations 
for many species. 



The model used here builds on its predecessors. Termed 
the "Gas/Surface" model it is a single point, gas and surface 
phase, chemical reaction network model with 279 species and 
2968 reactions. It includes the grain surface mechainsm used by 
Hasegawa et al. (1992) and Hasegawa & Herbst (119931) as well 
as the desorption mechanisms discussed by Willacy & Millar 
(119981 ). It is also fully "real-time-dependent". Instead of the dif- 
ferential equation approach taken by Shalabiea & Greenberg 
(1995) each separate time step has its own set of physical pa- 
rameters which are fed into the model from a storage file at the 
beginning of each time step. This approach has several advan- 
tages. Firstly it reduces computational strain on the model since 
there is no additional equation solving necessary. This is particu- 
larly important in situations where parameters are changing very 
quickly and there are major differences between adjacent time 
steps. (For example in the final stages of a cloud collapse where 
density increases drastically in a short time.) Secondly, this ap- 
proach allows the same set of software to be used to model mul- 
tiple different situations with no reprogramming at all. Only the 
input parameter file has to be changed. This allows great flexibil- 
ity in the scope of scenarios that could be investigated. Further, 
in certain situations where there are major, abrupt changes in 
physical parameters (eg. shocks), the differential equation ap- 
proach is highly prone to breakdown as generating a numerical 
solution for a given parameter requires some level of continuity 
between one time step and the next, wheras the parameter file 
approach avoids this problem completely. 



The Gas/Surface model is used here to investigate the sce- 
nario of complex molecule formation in grain mantles and their 
subsequent release back into the gas phase. We consider clouds 
at a variety of initial densities with different exposure to pho- 
toionization. As these clouds collapse to become denser and 
darker, their chemistry becomes more complex and their tem- 
perature drops allowing grain mantles to form. Simple gas-phase 
dominated chemistry at the beginning leads to the deposition of 
a variety of basic molecules onto an initially bare grain surface. 
Once on the grain, surface processes allow the build up of com- 
plex molecules in a frozen mantle. Later, as the cloud heats up, 
the mantle evaporates introducing the complex molecules into 
the gas phase with a variety of chemical consequences signifi- 
cantly different from gas phase only processing. 
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2. The Gas/Surface Model 

2.1. Design of the model 

The Gas/Surface model is used here to model a collapse sit- 
uation. A cloud collapse begins with the cloud optically thin 
and open to heating by external photons. Self shielding of H2 
is treated by the method of Wagenblast (1992 ). As it contracts 
it eventually becomes optically thick to external photons and its 
inner regions cool. It is this cooling which allows molecules to 
freeze out onto dust grain surfaces and form mantles. Eventually 
continuing collapse of the cloud, assisted possibly by the ignition 
of a proto-star at its center, causes the temperature to rise again 
desorping the mantles off the grain and returning mantle pro- 
cessed species to the gas phase. The model allows interchange 
between gas and surface phases, both freezing, thermal desorp- 
tion and a variety of other continously acting desorbtion mech- 
anisms, as described by Willacy (I1993I I and Willacy & Millar 
( |199 8i (See Section 1231 ) Once on the grain surface chemistry 
takes place as described in Section IZS] 

The reaction networks are solved numerically using the 
GEAR package (Hindmarsh 1 1 9741 [T983T l . It is the combination 
of both a time variable parameter space and interchange between 
gas and surface phase chemistry that makes the Gas/Surface 
model appropriate to investigate complex molecule formation. 
A fuller description of the model design can be found in Hall 
(fT997l l. 



Table 2. Mobile surface species 



Species 



Adsorption 
Energy Ep/k (K) 



Ref. 



1 


*H 


350.0 


1 


2 


*H2 


450.0 


1 


3 


*He 


100.0 


2 


4 


*C 


800.0 


1 


5 


*N 


800.0 


1 


6 


*o 


800.0 


1 


7 


*s 


1100.0 


1 


8 


*CH 


645.0 


3 


9 


*NH 


604.0 


3 


10 


*0H 


1259.0 


3 


11 


*CH2 


956.0 


3 


12 


*NH2 


856.0 


3 


13 


*CHi 


1158.0 


3 


14 




1107.0 


3 



References. (1) Tielen s & Allamandola ;i987i (2) T ielens & Hagen 
[T9821 (3) Allen & Robinson fTWTI 



where T is the grain surface temperature and Eh is the poten- 
tial energy barrier between adjacent surface sites, ^ Q.3Eo. 
The parameter vq is the characteristic vibration frequency for the 
mobile species (assumed isotropic). Its value (in Hz) is given by 



2.2. Surface chemistry mechanism 

The surface chemistry mechanism used in the Gas/Surface 
Model is taken from Hasegawa et al. (119921 ) and Hasegawa & 
Herbst (il993i with additional complex molecule formation reac- 
tions listed in the Appendix. Surface species are prefixed "*" to 
distinguish them from their gas phase counterparts. This scheme 
allows fourteen low molecular mass species (Table |2]i to be mo- 
bile on the grain surface while all other surface species are sta- 
tionary, held at "sites" on the surface. The adsorption energy Eo 
(in K) is the energy needed to liberate the species from the grain 
surface back into the gas phase. Parameter values (and their ref- 
erences) are listed in Table |4] 

The mobile species migrate around the grain surface by 
"hopping" from site to site. A mobile species arriving at an occu- 
pied site can react with either a stationary species already there 
or another mobile species which has also just arrived. Mobile 
species are regarded as occupying any given site for a finite pe- 
riod of time, generally about 10^'^s, the exact time being de- 
termined by species mass and prevailing temperature. There is 
an energy potential barrier separating adjacent sites and mo- 
bile species must possess sufficient energy to overcome it. For 
most of the mobile species "classical" type motions only are 
permitted. The species must have sufficient classical kinetic en- 
ergy to cross potential barriers. For atomic and molecular hydro- 
gen however quantum tunnelling through barriers is permitted. 
Which reactions take place is governed by a reaction set (Section 
I2.4I 1 analogous to those used for gas phase reactions, each with 
an accompanying rate coefficient. The grain itself is assumed to 
be inert and all surface chemistry involves the mantle only. The 
classical mobile species migrate over the grain surface. The time 
to move between two adjacent surface sites ("hopping time") ti^p 
is given by : 



vo = 



2n,Er 



(2) 



where «s is the surface density of sites (^ 1.5 x 10 m ) and m 
is the mass of the mobile species (in kg). This gives values of vq 
in the range lO'^-lO'^ s Hydrogen atoms and molecules are 
held to be quantum particles and their hopping times are given 
by the time for them to quantum tunnel through the potential 
barrier : 



1 

^hop 



= voexp 



—2a I 

— — y^lniEt 
n 



(3) 



10 



where a is the separation between two adjacent sites (^ 10 
m). The diffusion time, f,/,//, for a mobile species to scan the 
entire grain surface is given by : 



'diff 



Nsth 



nop 



(4) 



where A^, is the total number of surface sites on a grain (A^, = 
lO*"). The surface reaction rate coefficient kjj between two sur- 
face species / and j by classical diffusion is given by : 



kij - 



(5) 



where f,- and tj are the hopping times for species / and j and rig is 
the number density of grains. A gas:grain ratio of 100: 1 by mass, 
1:1.33 X 10"'^ by number is used. The factor k/j is the reaction 
probability for a specific reaction. /c,y =1.0 unless the associated 
activation energy (Ea) is non-zero. Non-zero activation energy 
enables certain reaction channels to be favoured over others. In 
the classical case : 



J_ ^ ^^g^ l-Eh\ 
thop \ kT I 



(1) 



exp 



(^) 



(6) 
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In the quantum case where at least one of the reactants is a 
hydrogen atom or molecule : 



Kij - exp 



—2a 



(7) 



where fi (kg) is the reduced mass in a two-body collision. 



2.3. Desorption mechanisms 

Desorption mechanisms can be divided into two basic classes, 
continuously acting which remove a few atoms/molecules at any 
one time and irregular such as shocks, which occur at random in- 
tervals but return most or all of a grain mantle to the gas phase. 
In the scenarios modelled here there are no shocks or similar 
mechanisms of sufficient strength to desorb mantles and con- 
tinous mechanisms dominate. Willacy (119931 1 provides an in- 
depth treatment of continuous mechanisms and four of them 
have been included in the models presented here : 

1 . Mantle Explosion 

2. Direct UV Photodesorption 

3. Cosmic Ray Induced Photodesorption 

4. Direct Heating by Cosmic Rays 

For each desorption mechanism an equation is provided in 
the reaction set with its own rate coefficient k : 

Surface species + desorption mechanism —> gas species (8) 

One of these four mechanisms, direct UV photodesorption, 
dominates when a cloud is in the early stages of collapse and 
still optically thin. Once collapse renders the bulk of the cloud 
opaque to external photons its effects are negligible and the 
other three mechanisms are more significant. However in the fi- 
nal stages of collapse direct thermal desorption of grain man- 
tles takes place and this outstrips all four continuous desorption 
mechanisms. The rate of evaporation is given by (Hasegawa et 
al. fT992b : 



= voexp 



(^) 



(9) 



For those species that are not mobile on the surface default val- 
ues are vo = lO'^s-' and = 2, OOOK. 



*CH30CH3. The reactions listed by Hasegawa et al. contain a 
number of laboratory measured potential energy barriers which 
slow down the rate of formation of particular species. Only a few 
such barriers are used here. Further, the reaction set includes al- 
most entirely reactions orientated to the formation of complex 
molecules. More reactions to produce other species would re- 
duce the complex molecule abundance by draining "raw mate- 
rial". It is recognized that this orientation of the reaction set is 
a limitation of the Gas/Surface model. However, it is regarded 
as a valid test of the hypothesis of complex molecule formation 
on grain surfaces. If, given all these advantages, such molecules 
still could not be produced by grain surface catalysis then the 
original supposition of such a production method would have to 
come under deeper scrutiny. 

A small number of surface reactions are included that are 
not directly involved in complex molecule formation. These are 
included because they are are believed to be significant to overall 
surface chemistry. The most basic of all is the grain formation of 
H2, whose equation is : 



*H + *H 



H2 



(10) 



It is considered that the H2 molecule produced passes di- 
rectly into the gas phase without the need for any specific desorp- 
tion process. This is consistent with quantum calculations which 
indicate the resultant molecule is formed with sufficient energy 
to eject it directly from the surface. Other reactions form small 
molecules which are frequently seen in high abundances in re- 
gions where grain mantles are thought to have recently evapo- 
rated. The basis of nitrogen chemistry is formation of *NH3 by 
the route : 



*N -H *H *NH H- *H ^ *NH2 + *H 



*NH3 



(11) 



Direct combination of atoms can also form *CN, *NO, 
*HNO, *HCN and *HNC. *02 and *OH can also form directly 
and *OH can further react to form *H20. Other combination re- 
actions allow the formation of *CO and *CH4. *CO and *H can 
also react to form *HCO and *H2 : 

*CO -H *H ^ *HCO -H *H ^ *CO+ *H2 (12) 

*H2CO can be formed and destroyed by the reactions ; 
*CH2 + *0 ^ *H2CO (13) 



2.4. Reaction set 

The basic gas phase reactions are chosen to give a representa- 
tive model of an interstellar cloud. Beginning with the UMIST 
astrochemical reaction rates database, RATE95 (Millar et al. 
11997 ), all reactions involving chlorine, phosphorus, iron, sili- 
con, sodium and sulphur are removed along with those involv- 
ing species with five or more carbon atoms. This yields a re- 
action set based on hydrogen, helium, carbon, oxygen and nitro- 
gen. Magnesium reactions are included to ensure a "token metal" 
and source of electrons. Additional provision is made for species 
formed only on grain surfaces which can pass into the gas phase. 
An extra set of reactions allowing gas phase destruction of these 
species has been added to the gas phase reaction set, listed in the 
Appendix (Table [TSTl. 

The surface phase reaction set ( Table [TSll is constructed with 
reference to those reactions provided by Hasegawa et al. ( 1992) 
and Hasegawa & Herbst ( 119931 ). However it is very much op- 
timized to allow the production of *CH30H, *C2H50H and 



*H2CO + *H ^ *HCO-H *H2 (14) 

There are two possible routes to *C02, which is stable on the 
grain surface and does not react further : 

*C0 + *OH *C02+ *H (15) 

*HCO 4- *0 ^ *C02-H *H (16) 

Direct combination of *CO with *0 has been omitted as it 
is believed that this reaction is extremely inefficient (Grim & 
d'Hendecourt 1986; Breukers |199I| l. Direct combination of car- 
bon atoms to form *C2 has also been omitted to simplify the sur- 
face chemistry. If it occurred it would lead by hydrogenation to 
the presence of significant amounts of *C2H6, not currently de- 
tected. C2H6 is seen in the cometry coma of comets Hyakutake 
and Hale-Bopp (Mumma et al. 1996) and it is not known if this 
formed on the surface or in the gas phase. The bulk of surface 
reactions comprise routes to the complex molecules *CH30H, 
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*C2H50H and *CH30CH3. *CH30H is formed by combining 
*CHn and *OH (n = 0-3) : 



*CHn + *OH ^ *CH„OH 



(17) 



The resultant is then hydrogenated as necessary to form 
*CH30H. Alternative routes are provided by : 



*CH3 + *0 



*CO + *H 



*CH30+ *H 



*COH 



*CH30H 



(18) 
(19) 



*COH can then hydrogenate. Reaction between *C0 and *H 
can also produce *HCO and the ratio of the two reactions is con- 
trolled by activation energy barriers. Routes to *C2H50H and 
*CH30CH3 involve forming a "spine" of C-C-O or C-O-C and 
then hydrogenating until the molecule saturates. This is accom- 
plished by the reactions : 



*CO -I- *C 



*CO + *C 



*CCO 



*COC 



(20) 
(21) 



Activation energy barriers can be included which slow down 
particular reactions. Some are from laboratory measurements. 
(See Hasegawa et al. ( |19921 i and references therein.) Those for 
the production of *C02 (eqns.[T5]& [16] above), *CCO (eqn.l20l 
above) and *COC (eqn.|2T| above) have been set to prevent ma- 
jor overabundances of these species. A full list of all surface 
reactions, with applicable activation energies, is given in the 
Appendix (Table [TSl l. It is recognized that the surface reaction 
set used is an approximation. It is constructed by interpolating 
between measured surface reactions, known or suspected pro- 
cesses in surface chemistry and the aimed-for result of producing 
complex molecule abundances and a grain mantle composition 
compatible with observations. 

2.5. Species set and starting abundances 

The species set is derived from the reaction set. (Full listing in 
the Appendix, Tables [16] and [T7]) All species can exist in the 
gas phase while only neutral species exist on the grain surface. 
Any ion freezing out is assumed to be neutralized. For ions with 
no neutral equivalent it is assumed they break up in exactly the 
same way as by electron neutralization in the gas phase. The 
starting abundances (Table [3]l are from Hasegawa et al. (!1992i 
Originally used by Leung et al. (119841) they are believed to be 
a reasonably accurate representation of the abundances found in 
an interstellar cloud, particularly one liable to collapse into a hot 
core or hot corino. Earlier gas phase only models derived using 
them are consistent with observations within the known accuracy 
of the chemistry. 



Table 3. Gas/Surface model - starting abundances 



Species 


Abundance w.r.t. (H + Ht) 


H 


1.00 


He 


1.40 X 10^' 


C 


7.30 X 10^'' 





1.76 X 10" 


N 


2.14 X 10^"^ 


Mg 


7.00 X 10-' 



3. Model Parameters 

3.1. Caiculating parameters 

As used here only four physical parameters are varied to model 
cloud collapse : density, visual extinction, gas temperature and 
grain temperature. All other parameters can be regarded as fixed 
and their values are given in Table ]4] The definitive factor is 
density, which controls the speed of collapse. In all cases here the 
cloud is assumed to be spherical and of uniform density. This is 
recognized as an approximation, possible effects of non-uniform 
density in star forming clouds are discussed by Shu et al. ( 19931. 
For the cloud to begin to contract under its own self-gravity its 
mass must exceed the Jeans mass (My) for its density. For the 
cloud densities modelled, the masses and corresponding Jeans 
mass are shown in Table [5] which also gives the relevant free- 
fall times (tff) and radii. 

The cloud collapse is halted at a density of 10^ cm"^, ob- 
servations of hot cores and hot corinos yield densities 10^-10^ 
cm"^. As the data in the Appendix demonstrate, by the end of 
the collapse all the hydrogen in the cloud has become molecular. 
(H2 abundance = 5.0 xlO"' with respect to total hydrogen.) Ay 
is calculated from Brown ( 1988) and its increase is stopped at 60 
magnitudes as further increases in have negligible effect on the 
chemistry. 

All cloud collapse models presented here assume that gas 
and grain temperature are the same. Since the gas and grains are 
closely coupled and the temperature mainly varies as a smooth 
function it is unlikely that they would be significantly differ- 
ent over the temperature region where the grains have mantles. 
As a cloud collapses gravitational energy liberated can be effi- 
ciently radiated in the infrared and submillimeter bands (Dyson 
& Williams 1980). However, at some point the temperature must 
rise to begin nuclear processing. It is assumed that below a den- 
sity of 10^ cm""* the cloud can efficiently radiate and temperature 
(T) in this region is found from a curve fitting technique derived 
by Tarafdar et al. (iT985T l : 
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2.5 + log An) - logeil + 500ex/?(-1.8Ay)) 



(22) 



where n is number density and Ay visual extinction. Above 
n = 10^ cm the temperature is determined by the release of 
gravitational energy : 



(23) 



where Ta is the temperature at 10^ cm"-'. From Tarafdar et al. 
To - 10 K and no = 10^ cm""*. The final temperature obtained is 
T - 46.4 K at the maximum number density n = 10^ cm"-'. The 
final environments produced by the models have a temperature 
of r = 46.4 K, a number density n = 10^ cm"^, an Ay of 60 mags 
and all remaining hydrogen molecular. 

Recently a number of "hot corinos" (Ceccarelli et al. 120041 
Bottinelli et al. I2004I I have been investigated. These are re- 
gions of low mass (^1 Mq) star formation embedded within 
and condensed from molecular clouds. They have typical sizes 
of <200 AU, temperatures of ^100 K and densities of >10^ 
cm"^ (Cazaux et al. 2003, Bottinelli et al. 2004 and refer- 
ences within.) Complex molecules have been detected in four 
hot corino sources, listed in Table [T] 
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Table 4. Gas/Surface model - 
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fixed parameters 



Parameter 


Value 


Ref. 


UV Radiation field 


3.0x10** cm"-s"'Hz"' 


1 


Cosmic ray ionization rate 


1.3xlO""s"' 


2 


Grain albedo 


0.5 


3 


Sticking coefficient 


1.0 


4 


Sticking coefficient for hydrogen atoms 


U.i 


c 

J 


Yield per UV photon on grain mantle impact 


IQ-* 


6 


Surface density of sites 


1.5x 10'^ cm-2 


4 


Number of sites per grain 


10* 


4 


Fractional density of grains (w.r.t. Density of H) 


1.33x10^'^ 


4 


Average grain radius 


10"^ cm 


4 


Barrier width between adjacent sites 


10-** cm 


4 


Grain velocity (non-thermal) 


10" cms-' 


7 



References. (1) Draine il978i (2) van Dishoeck & Black a986i (3) Whittet 1992 , (4) Hasegawa et al. il992l (5) Brown [T988] (6) Willacv [T993] (7) 

Volk et al .l980. 

Table 5. Gas/Surface model - cloud parameters 



Initial Number 
Density (cm^^) 


Mass 
(Me) 


Mj 

(Mo) 


iff 

(years) 


Initial Radius 
(parsecs) 


Final Radius 
(parsecs) 


Initial 
A„ (Mag) 


Initial 
Temp. (K) 


100 


1700 


1692.7 


4.40x10" 


5.0818 


0.065 


1.39 


48.5 


1,000 


150 


144.3 


1.39x10" 


1.0501 


0.034 


2.86 


20.2 


10,000 


26 


25.97 


4.40x10^ 


0.2717 


0.006 


7.41 


13.9 



Table 6. Collapse scenario model parameters 





Initial Number 


Parameter 


Minimum 


Plateau 




Density (cm^^) 




Temperature 


Point 


Model 1 


100 


Time (yrs) 


4.36 X 10" 


4.39 X 10" 


Free-fall 


Atomic 


n (cm^^^) 


6.97 X 10"* 


1.00 X 10** 


Model 2 


1,000 


Time (yrs) 


1.36 X 10" 


1.39 X 10" 


Free-fall 


Atomic 


n (cm^^^) 


6.78 X 10' 


1.00 X 10** 


Model 3 


10,000 


Time (yrs) 


4.18 X 10" 


4.38 X 10' 


Free-fall 


Atomic 


n (cm"^^) 


9.47 X 10*^ 


1.00 X 10** 


Model 4 


100 


Time (yrs) 


2.18 X 10' 


2.19 X 10' 


Retarded (x 5) 


Atomic 


n (cm^^) 


1.00 X 10" 


1.00 X 10** 


Model 5 


1,000 


Time (yrs) 


6.84 X 10" 


6.93 X 10" 


Retarded (x 5) 


Atomic 


n (cm^^^) 


9.87 X 10" 


1.00 X 10** 


Model 6 


10,000 


Time (yrs) 


2.07 X 10" 


2.19 X 10" 


Retarded (x 5) 


Atomic 


«(cm"'') 


7.73 X lO"* 


1.00 X 10** 


Model 7 


10,000 


Time (yrs) 


4.18 X 10' 


4.38 X 10' 


Free-fall 


Molecular 


n (cm^^^) 


9.47 X 10' 


1.00 X 10** 


Model 8 


10,000 


Time (yrs) 


2.07 X 10" 


2.19 X 10" 


Retarded (x 5) 


Molecular 


n(cm"^) 


7.73 X 10' 


1.00 X 10** 



3.2. Collapse scenarios 



In all, eight different cloud collapse scenarios are investigated, 
designated models 1-8. Models 1-3 are free-fall collapses with 
starting abundances of 100, 1,000 and 10,000 cm"^, regarded as 
slow, medium and fast. The observed formation rate of new stars 
is too low for all stars in the galaxy to be formed by free-fall 
collapse. It is widely asserted that some mechanism supports 
clouds against collapse, the most commonly suggested being 
magnetic fields (Shu et al. 1993 ). For this reason models 4-6 are 
retarded collapses, again with starting abundances of 100, 1,000 
and 10,000 cm"^. No specific retardation mechanism is incorpo- 
rated in the Gas/Surface model nor is one necessary. Retardation 
is accounted for by arbitrarily slowing down the decrease in ra- 
dius by a factor of 5.0. This ensures that the cloud takes longer to 
collapse with corresponding further evolution in the chemistry. 
Apart from retardation the parameter change in density, gas and 



grain temperature and visual extinction follows exactly the same 
pattern as in models 1-3. 

For the purposes of direct comparison as many parameters 
as possible are kept constant for all collapses. At a density of 
100 cm"^ the hydrogen is likely to be almost completely atomic. 
Observed clouds of this density are optically thin and exposure 
to background starlight ensures an ample supply of UV photons 
to break up H2. At 1,000 cm"^ this is still approximately true 
though some H2 does form. However at the high density case of 
10,000 cm"^ it is unphysical. The clouds are sufficiently dense 
that they will be predominantly H2. For this reason the final mod- 
els, 7 and 8, replicate the free-fall and retarded collapses from 
10,000 cm with their initial hydrogen abundance molecular 
instead of atomic. 

Essensially the chemistry during cloud collapse can be di- 
vided into three approximate segments : 
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1. Pre-freeze out. Gas phase chemistry dominates, similar to 
that in gas phase only models. 

2. Freeze out. As the temperature falls and the visual extinc- 
tion (Ay) increases (curtailing photodesorption), freeze out 
occurs and surface reactions dominate the chemistry. This 
occurs around the minimum temperature, ^ 10.0 and 10.2 K. 

3. Post- freeze out. Once the temperature rises the entire grain 
mantle passes back into the gas phase. Grain surface prod- 
ucts are suddenly released into a very dense, warm, dark 
cloud environment. This leads to a rapidly changing chem- 
istry where a wide variety of different species interact. 

Table |6] lists the parameter values for each model at the be- 
ginning of its run, the point of minimum temperature and the 
"plateau point". The plateau point is where parameter change 
ceases. In all cases at the plateau point the density is 10** cm"^^, 
temperature 46.4 K and visual extinction 60.0 magnitudes. For 
the 2.0 X 10^ years beginning at the plateau point each model 
most closely represents the chemistry of a hot corino. All model 
runs terminate 2.0 x 10^ years after the plateau point is reached. 

4. Discussion 

The Appendix contains output data from each of the eight mod- 
els. The most abundant gas phase species at both the plateau 
point and 2.0 xlO^ years later are listed in Tables [19] - l26l 
Variation in gas phase species abundance over time is shown in 
Figures ?? -[8] The results are discussed in this section. 

4. 1 . Mantle composition 

The test of any theoretical model is to compare it to observations. 
Table |7]presents observational data on interstellar ices towards a 
number of sources. The abundances shown are relative to *H20 
and have been measured by Boogert et al. ( 2008 ) for *HCOOH 
& *CH30H, Pontoppidan et al. (2008} for CO & CO2, Oberg et 
al. (|20D^ for CH4 and Gibb et al. (2000^ for *NH3 & *H2CO. 
For data taken from Boogert et al. (2008 ) measured values have 
been taken, observations where a lower limit only has been es- 
tablished have been ommitted. Further Boogert et al. (.2008i) data 
for *NH^ has been ommitted as the authors themselves point out 
that the inclusion of certain components in the abundance esti- 
mation is still a matter of debate. 

It is assumed that ice measurements reflect the composition 
of grain mantles along the measured lines of sight. No measure- 
ments are available for mantles in hot cores or hot corinos as 
within these regions the temperature is sufficiently high that any 
grains have been fully thermally desorbed. Observations made 
of other, cooler regions are therefore the best available indicator 
of grain mantle composition. All models produce mantles with a 
predominance of *H20, (absolute composition 40-64%), agree- 
ing with observations. Table [8] presents mantle composition at 
greatest thickness as a percentage of *H20 abundance for all 
mantle species with abundance > 1% *H20 and for the com- 
plex molecules regardless of their contribution. Table |9] shows 
the conditions during each model when the mantle is at its great- 
est thickness. The point of greatest thickness is defined as that 
timestep with the largest number of surface layers. Where this is 
the same for two or more timesteps the latest is taken. 

*CO is well modelled in models 1-4 and 6 and overproduced 
in the other three. For models 7 and 8 this is a consequence of 
the molecular hydrogen starting abundance. In both cases atomic 
carbon and oxygen are able to preferentially combine with each 



other due to the absence of atomic hydrogen which would oth- 
erwise convert them to *H20 and *CH4. Model 5 appears to be 
an optimum case (with initial atomic hydrogen) for the produc- 
tion of *CO and also subsequently *C2H50H and *CH30CH3, 
hence the higher proportion of all three of these species. Table 
[TO]bears this out, showing an order of magnitude lower carbon 
abundance at freeze out time in model 5 compared to the other 
models. Correspondingly both *CO and *C02 are significantly 
more abundant in model 5. Both these are preferentially formed 
in the gas phase before freeze out. *C02 is well modelled in 
models 4, 5, 7 and 8 and underproduced in all others. Again in 
models 7 and 8 the absence of initial atomic hydrogen is the ma- 
jor factor here, while in models 4 and 5 it seems that the longer 
time period before freeze out enhances *C02 production, for the 
same reasons as it also enhances *CO production. Exactly the 
same factors result in *CH4 being well modelled in models 5, 7 
and 8 and overproduced in all others. The faster collapses (mod- 
els 1, 2, 3 and 6) freeze out atomic carbon before it can form 
molecules in the gas phase and once on the grain this rapidly 
hydrogenates into *CH4. Model 4 then appears as a transitory 
case. It has more *CH4 than model 5, about the same as the 
faster cases, with correspondingly more *CO and *C02 than the 
faster cases and less than model 5. It appears that for model 4 
the greater time before freeze out for CO and CO2 to form is 
partially offset by the lower starting density. The same effect can 
be seen when comparing models 1 and 2, though it is not as pro- 
nounced. 

*NH3 is well modelled in models 1, 2 and 3, above the 
1 % threshold in model 7 and underproduced in the other mod- 
els. Here it seems that rapid collapse alone is the major factor 
Abundance in model 7 is significantly lower than in models 1-3, 
entirely a consequence of the molecular hydrogen starting abun- 
dance. Models 4-8 all predict > 1% *HNO. Correspondingly 
these are the models with the least *NH3. The highest *HNO 
abundance is in models 7 and 8, again a consequence of the 
molecular hydrogen starting abundance. Nitrogen cannot easily 
form NH3 in the gas or surface phase due to the lack of atomic 
hydrogen. Instead, atomic nitrogen freezes out and combines 
with *OH to form *HNO. The lack of atomic hydrogen allows 
there to be proportionately more *OH available, which other- 
wise would hydrogenate to form *H20. Models 4 and 5 seem to 
be intermediate cases, their longer collapse times allowing some 
but not all of the initial gas phase atomic hydrogen to convert to 
molecular before freeze out. 

*HNO is not currently observed in grain mantles and is an- 
other model prediction to be searched for. However its greatest 
abundance occurs in two collapses (models 7 and 8), both with 
high starting density. A dense cloud with ample molecular hy- 
drogen could have evolved chemically before collapse begins 
reducing the *HNO abundance by enabling atomic nitrogen to 
combine into other species in the early gas phase. All models 
predict a significant mantle abundance of *N2 and most models 
*02, at ^ 10%. In the absence of observational measurements 
these must be regarded as model predictions. Ehrenfreund & van 
Dishoeck d 19981 1 suggest possible ways of searching for *N2 and 
*02 while Ehrenfreund & Schutte (I2000I I describe the current 
state of observations. 

Although mantle composition observations cannot be re- 
garded as complete and definitive, approximate agreement for 
*H20 and *CO is seen in mantles produced by the Gas/Surface 
model. Mantle compositions in models 3 and 6 (10,000 cm"^ 
free fall and retarded) must be regarded as suspect as it is likely 
that the scenarios themselves are unphysical, starting as they do 
from atomic hydrogen. Models 4, 5, 7 and 8 produce composi- 
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Table 7. Observed interstellar ices (See text for references) 



Source 


*H20 


*HCOOH 


H-CHjOH 


*CH4 


*C02 


*C0 


*NH3 


*H2CO 


Bl-b 


100.00 


3.1 


11.2 


3.3 










IRAS 08242-5050 


100.00 


2.7 


5.5 


5.0 










IRAS 15398-3359 


100.00 


1.9 


10.3 


6.0 


35.24 


5.63 






SVS 4-5 


100.00 


4.6 


25.2 


6.1 


30.46 


13.09 






L1014 IRS 


100.00 


5.4 


3.1 


7.1 










W33A 


100.00 


5.2 


14.7 




12.49 


7.37 


15.0 


6.0 


GL2136 


100.00 


5.1 


8.5 




13.24 


10.19 






GL 7009S 


100.00 


2.5 


31.3 












NGC7538 IRS9 


100.00 




7.5 








13.0 


4.0 



Table 8. Mantle composition w.rt H2O at maximum mantle thickness for each model 



Species 


Model 1 


Model 2 


Model 3 


Model 4 


Model 5 


Model 6 


Model 7 


Model 8 


*H20 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


*CH4 


34.23 


28.62 


32.50 


28.34 




31.51 


3.08 


2.45 


*NH3 


5.94 


3.88 


5.04 








1.14 




*N2 


3.75 


5.60 


4.39 


7.87 


13.76 


6.31 


8.46 


10.50 


*CO 


3.65 


7.15 


3.88 


14.73 


69.53 


4.42 


55.66 


71.91 


*C02 




1.69 




7.16 


24.46 


2.22 


22.12 


28.43 


*02 


1.24 






5.65 


12.79 


2.17 


2.77 


15.49 


*HNO 








1.92 


3.71 


1.45 


9.16 


11.68 


*C4H 
















1.17 


*CH30H 


7.98 


14.70 


11.46 


10.32 


12.12 


13.05 


1.91 


1.57 


*C2HjOH 


1.19 X 10^' 


7.74 X 10^* 


7.64 X lO^"* 


9.50 X 10^* 


1.28 X 10^' 


6.21 X 10" 


3.49 X 10"- 


2.43 X 10^^ 


*CH30CH3 


5.80 X 10-' 


2.60 X 10^' 


2.10 X 10^' 


9.00 X 10^** 


9.47 X 10^^ 


5.50 X 10-' 


2.24 X 10^^ 


1.05 X lO-'* 



Table 9. Mantle layering 







Model 1 


Model 2 


Model 3 


Model 4 


Maximum 


Layers 


184 


169 


179 


154 


Thickness 


Time (years) 


4.38 X 10^ 


1.38 X 10*^ 


4.28 X 10' 


2.19 X 10' 




Density (cm^^) 


5.95 X 10^ 


5.02 X 10*^ 


3.39 X 10^ 


5.62 X 10*^ 




Temperature (K) 


18.1 


17.1 


15.0 


17.8 






Model 5 


Model 6 


Model 7 


Model 8 


Maximum 


Layers 


121 


169 


126 


119 


Thickness 


Time (years) 


6.90 X 10** 


2.14 X 10*^ 


4.28 X 10' 


2.14 X 10* 




Density (cm^^) 


8.81 X 10^ 


3.20 X lO'' 


3.39 X 10« 


3.20 X 10*^ 




Temperature (K) 


20.7 


14.7 


15.0 


14.7 



tions closest to observations, although models 7 and 8 predict 
*CH30H abundance at the low end of its observed abundance 
range. The effects of collapse timescale on mantle composition 
are shown in Table [TO] which lists relative abundances of gas 
phase C, CO, CO2, H, H2 and O with respect to total hydrogen 
at the beginning of freeze out. This is arbitrarily defined, for the 
purposes of this discussion, as the last timestep before at least 
five mantle layers form on the grain. Models 2, 4 and particu- 
larly 5 seem to have an optimum combination of collapse speed 
and density to allow significant CO to form in the gas phase. 
This partially depletes atomic carbon and enables freeze out to 
occur while there is still some atomic hydrogen available with 
coiTesponding consequences for the chemistry thereafter. 

4.2. Complex molecules 

4.2.1. Complex molecules - surface phase 

The prime function of the Gas/Surface model is to investigate the 
formation of complex molecules by grain surface catalysis. Gas 
phase abundances in the final stages of collapse are governed 
by surface abundances prior to mantle desorption. The mantle 
abundances with respect to *H20 of *CH30H, *C2H50H and 



*CH30CH3 at the greatest mantle thickness for each model are 
Usted in Table[8] Models 1-6 all produce *CH30H in the 5-15% 
range, agreeing well with observations. Models 7 and 8 pro- 
duce ^1%. This is a consequence of the initial hydrogen being 
molecular *CH30H requires atomic hydrogen for its formation 
(Section [2.41 ) and with little available formation is significantly 
reduced. 

Models 1-6 all produce about the same amount of *C2H50H 
(^ 0.001%) while models 7 and 8 produce significantly more 
(^ 0.025%). For *CH30CH3 models 1-4 and 6 produce ^ 
0.0000002% while models 5, 7 and 8 production is higher at 
- 0.0001%. Formation of complex molecules requires atomic 
carbon and *CO on the grain, the *CO originating in the gas 
phase (Section 12741 1. Maximum production takes place when the 
prior gas phase has produced significant CO, while leaving suf- 
ficient atomic carbon to form C-C-O and C-O-C structures on 
the surface after freeze out. At the same time if there is consider- 
able atomic hydrogen present at freeze out then this will rapidly 
hydrogenate atomic carbon on the grain surface leaving little 
available to form complex molecules. This is demonstrated by 
models 7 and 8 where the initial hydrogen is entirely molecular. 
Absence of much atomic hydrogen leaves CO to form preferen- 
tially in the gas phase. It then freezes out alongside atomic car- 
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Table 10. Gas phase abundances - freeze out 
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Time (years) 


C 




CO 




CO2 




H 


H2 







Model 1 


3.67 xlO" 


6.76 xlO" 





8.34 xlO" 


^7 — 


7.95 xlO" 


~n — 


7.79 xlO"' 


1.11 xlO"' 


1.71 xlO 


— 


Model 2 


7.94 xlO*^ 


5.81 xlO" 


-5 


1.27 xlQ- 


-5 


3.97 XlO- 


-11 


7.64 XlO"' 


1.18 XlO"' 


1.61 XlO 


-4 


Model 3 


1.25 xlO^ 


6.99 xlO- 


-5 


9.23 xlO 


-7 


1.07 XlO- 


-12 


9.73 XlO"' 


1.33 xlO"- 


1.71 XlO 


-4 


Model 4 


1.81 xlO' 


5.67 xlO" 


-5 


1.13 xlO 


-5 


1.46 XlO 


-11 


3.10 XlO"' 


3.45 XlO"' 


1.60 XlO 


-4 


Model 5 


4.42 xlO^ 


3.69 xlO- 


-6 


6.55 xlQ- 


-5 


3.69 xlO- 


-7 


2.20 XlO"' 


3.90x10-' 


9.64 XlO 


-5 


Model 6 


1.25 xlO"* 


6.99 xlO- 


-5 


9.23 xlQ- 


-7 


1.07 xlO- 


-12 


9.73 XlO"' 


1.33 xlQ-^ 


1.71 XlO 


-4 


Model 7 


2.51 xlO'* 


6.08 xlO- 


-5 


6.07 XlO- 


-6 


5.95 xlO- 


9 


3.48 xlO"^ 


5.00x10-' 


1.61 XlO 


-4 


Model 8 


2.51 xlO"* 


6.08 xlO- 


-5 


6.06 XlO- 


-6 


5.94 XlO- 


9 


3.47 XlO"' 


5.00 XlO-' 


1.61 XlO 


-4 



Table 11. Complex molecule abundances - gas phase 











Abundance w.r.t. H + H2 












CH3OH 


C2H5OH 




CH3OCH3 




Plateau Point 


Plateau Point 


Plateau Point 


Plateau Point 


Plateau Point 


Plateau Point 








+ 2.0 xlO^ yrs 




+ 2.0 XlO' yrs 






+ 2.0 XlO' yrs 


Model 1 


1.12 xlO- 


-5 


4.32 XlO-' 


1.78 XlO-" 


1.45 XlO-** 


3.01 XlO 


-7 


8.36x10-' 


Model 2 


1.83 xlO- 


-5 


3.80x10-' 


1.04 XlO-' 


3.32 XlO-** 


4.11 xlO- 


-7 


5.81 XlO-' 


Model 3 


1.55 xlO- 


-5 


5.23 XlO-' 


1.11 XlO-' 


1.26 XlO-** 


3.63 xlO- 


-7 


1.04 xlO-« 


Model 4 


8.56 xlO- 


-6 


7.98 XlO-** 


1.69 XlO-' 


6.21 XlO-** 


8.61 xlO- 


-7 


6.86 XlO-** 


Model 5 


5.44 xlO- 


-6 


1.34x10-' 


1.40 XlO-' 


6.48 XlO-'" 


7.69 xlO- 


-7 


1.62 XlO-" 


Model 6 


1.09 xlO- 


-5 


1.66 XlO-' 


1.76 XlO-' 


5.63 XlO-' 


1.64x10- 


-6 


1.92 XlO-' 


Model 7 


8.71 xlO- 


-7 


2.72x10-' 


2.30 XlO-** 


4.76 XlO-'" 


1.94x10 


-7 


7.07 XlO-" 


Model 8 


1.51 XlO 


-7 


1.97 XlO-' 


9.32 XlO-' 


2.73 XlO-'" 


1.27 xlO- 


-7 


3.94 XlO-" 



Table 12. Molecular abundance in a sample of four low-mass protostars (BottinelU et al. l2007l l 



Molecule 


IRAS 16293-2422 


NGC1333-IRAS4A 


NGC1333-IRAS4B 


NGC1333-IRAS2A 


Ref. 


H2CO 


1.0 X 10-' 


2.0 X 10-** 


3.0 X 10-" 


2.0 X 10-' 


1,2,3 


CH3OH 


1.0 X 10-' 


< 1.0 X 10 ** 


7.0 X 10-' 


3.0 X 10-' 


4 


HCOOCH3-A 


1.7 X 10-' 


3.4 X 10"* 


1.1 X 10-** 


< 6.7 X 10-' 


5,2,6 


HCOOH 


6.2 X 10-** 


4.6 X 10-' 


< l.Ox 10-" 


< 1.2 X 10' 


5,2,6 


CH3OCH3 


2.4 X 10-' 


< 2.8 X 10-** 


< 1.2x 10-" 


3.0 X 10-** 


5,2,6,7 


CH3CN 


1.0 X 10-** 


1.6 X 10-' 


9.5 X 10-** 


8.7 X 10-' 


5,2,6 


C2H5CN 


1.2 X 10-** 


< 1.2 X 10-' 


< 7.5 X 10-' 


< l.Ox 10-' 


5,2,6 



References. (1) Ceccarelli et al. 2000. (2) Bottinelli et al. 2004) (3) Maret et al. 2004; (4) Maret et al.l2005r 



(5) Cazaux et al. 2003, (6) Bottinelli et al. l2007l (7) J0rgensen et al. . 20051 



bon which cannot readily hydrogenate into *CH4 and so is avail- 
able to form C-C-O and C-O-C. Models 7 and 8 show low *CH4, 
enhanced *CO and the largest proportions of both *C2H50H and 
*CH30CH3. 

Models 1, 2, 3 and 6 with their faster collapse times and 
initial atomic hydrogen all have significantly lower *CO in the 
mantle and hence much lower complex molecule abundance. 
Models 4 and 5 have the longest collapse time of any of the mod- 
els. This leads to significantly more *CO, in the case of model 
5 approximately the same as in the molecular hydrogen start- 
ing abundance cases. Model 4 again appears as a transitory case. 
It has 15% *CO, formed in the gas phase and then frozen out. 
However its initial low density means that significant atomic hy- 
drogen is present at freeze out and this can hydrogenate the re- 
maining atomic carbon to *CH4 (28%) leaving very little for the 
production of complex molecules. Model 5 avoids this situation. 
Here substantial *CO is frozen out but very little atomic hydro- 
gen remains, hence the paucity of *CH4. Much of the remaining 
carbon has contributed to the high *C02 abundance (also seen 
in models 7 and 8) however some has led to the significantly 
increased *CH30CH3. 

Each model produces a higher mantle abundance for 
*C2H50H than *CH30CH3. For models 1-4 and 6 the differ- 
ence is greater than 3 orders of magnitude while for models 7 



and 8 it is 2 and for model 5 only 1 . The abundance difference 
seems to be caused by reaction channel differences. *C2H50H 
forms mainly by hydrogenation on the grain of C-C-O. However 
*CH30CH3 seems to form mainly by freeze out from the gas 
phase. Its gas phase origin is a two stage process with C2H60H^ 
(protonated dimethyl ether) as an intermediate ; 

CHJ -I- CH3OH ^ C2H6OH+ + photon (24) 

C2H6OH+ + electron CH3OCH3 + H (25) 

Gas phase C2H60H^ can also be destroyed by : 

C2H6OH+ + electron CH3 + CH3OH (26) 

Essentially C2H60H^ is made from CH3OH. Once formed 
it is then destroyed in forming CH3OCH3 (eqn. l25l l which can 
freeze out onto the grain surface where it is stable or it can 
form CH3OH (eqn. l26T l which is then available to form more 
C2H60H^ via eqn.|24] As long as the temperature is low enough 
to ensure that there is a net freeze out of molecules onto the 
grain surface this process results in a build up of *CH30CH3. 
While workable this gas phase production of CH3OCH3 which 
then freezes out is less efficient than the direct grain surface hy- 
drogenation which produces *C2H50H in significantly greater 
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abundances. Although production of CaHeOH^ takes place in 
the gas phase grain surface catalysis is crucial to the process. At 
the time of freeze out most CH3OH is being produced by evapo- 
ration from the grain mantle. In turn the dominant production of 
♦CHsOH is by a chain of direct hydrogenation reactions of *C0. 

Comparison of reaction rates of complex molecule formation 
for the different models shows that model 5 seems to have an 
optimum combination of parameters to form *CH30CH3, hence 
its significantly higher abundance in this model. Table [1] does 
not show any consistency in C2H5OH : CH3OCH3 ratio and it 
seems from the models here that the abundance ratio of these 
two is significantly sensitive to prevailing conditions. 

4.2.2. Complex molecules - gas phase 



Table 13. Acceptable gas phase fractional abundance ranges for 
complex molecules 



Molecule 


Minimum 


Maximum 




Abundance 


Abundance 


CH3OH 


5.0 X 10-** 


2.0 X 10-" 


CH3OCH3 


8.0 X 10-" 


7.4 X 10-' 



Other than *CH30H no observations are currently avail- 
able of complex molecules in grain mantles. However measure- 
ments do exist of gas phase abundances (Table [T2l i and these 
can be compared to the results of the models shown in Table 
nn Comparisons are made at the plateau point and at 2.0 x 10^ 
years later, the end point of time segment 2. Between these two 
points the model most closely resembles a hot corino. There are 
uncertainties in observations and even when these can be min- 
imized apparently similar sources yield significantly different 
abundance values for the same species. There are further uncer- 
tainties in the chemistry both in the processes of surface chem- 
istry and the accuracy of gas phase rate coefficients. For these 
reasons a model-predicted gas phase abundance is considered to 
be in reasonable agreement with observations if it falls within 
half an order of magnitude of the range of observed values at 
any time between the plateau point and 2.0 xlO^ years later. The 
acceptable ranges for complex molecules are given in Table [T3l 
(derived from Table [T2li. 

CH3OH is well modelled in all scenarios, though the abun- 
dance varies significantly. In all of models 1-6 (all models 
starting with atomic hydrogen) CH3OH is overproduced at the 
plateau point and then for five of these cases falls to within the 
acceptable range by 2.0 x 10^ years later. The sole exception here 
is model 5 (retarded collapses from 1,000 cm-^^) where the abun- 
dance actually drops to below the valid range 2.0 x 10^ years 
later Since the abundance has passed through the valid range 
during the time period in question this is held to indicate valid- 
ity. Models 7 & 8 (with molecular hydrogen starting abundance) 
are both within the valid range at the plateau point and fall be- 
low it by 2.0 X 10^ years later. In all cases the CH3OH gas phase 
abundance declines after the grain mantles are evaporated. 

CH3OCH3 can also be considered well modelled in all cases, 
though its behaviour is more diverse. Models 1, 2 & 3 are within 
acceptable range at the plateau point and by 2.0 x 10^ years later 
abundance has actually increased, suggesting further production 
in the gas phase. At this time model 2 is still within the accept- 
able range while 1 & 3 overproduce. Models 4, 5, & 6 all over- 
produce slightly at the plateau point while 2.0 x 10^ years later 
abundances have dropped to within acceptable range for models 



4 & 6 and below acceptable range for model 5. As with CH3OH 
model 5 CH3OCH3 abundance is considered to be validly mod- 
elled as it has passed through the acceptable range during the 
relevant time period. In this case the preciptious drop, 4 orders 
of magnitude, serves to illustrate how volatile the chemistry be- 
comes in the period after the grain mantles evaporate. Models 7 
& 8 show a similar behaviour being within acceptable range at 
the plateau point and also dropping by 4 orders of magnitude 2.0 
X 10^ years later 

As discussed in Sections 14.11 and 14.2.11 the production of 
complex molecules is dependent on the gas phase abundances 
of atomic carbon and hydrogen at the point of freeze out. Gas 
phase processes are considerably better understood than surface 
processes, and in particular the potential barriers used could be 
different leading to changes in the relative proportions of species. 
With these caveats we conclude that complex molecules, seen 
in both hot cores and hot corinos, can be made in the ob- 
served abundances by grain surface catalysis processes. Ethanol 
(C2H5OH) production has been included in the models here as 
it is frequently seen in the same environments as dimethyl ether 
(CH3OCH3) (Table[TJ. Currently it has not been observed in hot 
corinos and the abundances modelled here must be regarded as 
predictions awaiting the test of observation. The only major dis- 
tinction that can be draw between the different models is be- 
tween those starting with molecular hydrogen (models 7 & 8) 
compared to the other six. As would be expected the lack of 
initial atomic hydrogen stifles the production of those molecu- 
lar species which require it. Correspondingly this then allows 
higher abundances of other species whose "raw material" would 
otherwise have been removed by combination with atomic hy- 
drogen. No correlation is noticable in the production of complex 
molecules in either gas or surface phase with the cloud starting 
mass and/or collapse time. This may be indicative of a volatile 
chemistry and/or of the limitations of our current understanding. 

Table [12] lists observational measurements for a number of 
other molecules observed in hot corinos. These too can be com- 
pared to the models here as a validity check. HCOOH and 
C2H5CN are not in the Gas/Surface model reaction set as used 
here. 

Both of these two species are obvious candidates for inclu- 
sion in future modelling. HCOOH may be particularly signifi- 
cant as both gas (Table \T% and solid phase (Table |7]) data are 
now available for it. Both Bottinelli et al. (120071) and Fuente 
et al. (2005) point out that HCOOH shows significantly differ- 
ent abundance behaviour between corino type objects, as mod- 
elled here, and higher temperature hot cores, unlike a majority 
of other observed species. The reasons for this are not known 
but present an interesting issue for prospective future models. 
H2CO, HCOOCH3 and CH3CN are included and their abun- 
dances are listed in Table [15] The same criteria for acceptable 
range is used as above, namely a gas phase abundance is consid- 
ered to be in reasonable agreement with observations if it falls 
within half an order of magnitude of the range of observed val- 
ues at any time between the plateau point and 2.0 xlO^ years 
later Relevant values are shown in Table [T4l 

H2CO is well modelled in all collapse scenarios. Table [12] 
lists only abundance data for HCOOCH3-A. Bottinelli et al. 
(2007), from where this data is taken, state that the abundance 
of the "E" form of HCOOCH3 is usually very close to that of the 
"A" form and that they regard the HCOOCH3 total abundance 
to be twice that of HCOOCH3-A. For this reason the accept- 
able range of HCOOCH3 is given in Table[Tlas 1.8 xlO - 7.2 
xlO-^, i.e., half an order of magnitude above and below double 
the HCOOCH3-A abundance. Modelled HCOOCH3 is under- 



p. Hall and T. J. Millar: Complex molecule formation in grain mantles 



11 



abundant at the plateau point in all cases. By 2.0x10^ years later 
it remains underabundant in seven of the eight cases. Only in 
model 2 does it rise above the acceptable lower limit. CH3CN is 
well modelled in most of the cases. In models 1-4 & 6 it is within 
the acceptable range both at the plateau point and 2.0x10^ years 
later. In models 5 & 8 it is always underabundant. In model 7 it 
is within acceptable range at the plateau point and has fallen to 
below it by 2.0x10^ years later. For all three of these molecules 
there is no descernible correlation between abundance produced 
and cloud starting mass and/or collapse time. 



University Department of Applied and Engineering Physics for helpful conver- 
sations about optical absorbtion. We both wish to thank an anonymous referee 
whose comments improved an earlier version of this paper 



Appendix - Model Data 



5. Conclusion 

Overall the Gas/Surface model as used here reproduces to a 
reasonable degree the observed complex molecule abundances 
seen in hot corinos. In particular abundance values produced for 
methanol and dimethyl ether are comparable to those observed 
in both gas and surface phases. 

The surface phase molecular abundances produced by the 
model reflect the starting conditions, particularly the nature of 
the initial hydrogen - atomic or molecular This determines the 
degree of processing in the initial gas phase and hence the abun- 
dance of species frozen on to the dust and available for surface 
chemistry. This in turn aUows an understanding of the differ- 
ences in grain mantle composition observed in different regions. 
The model results clearly indicate that both the starting condi- 
tions and particularly the timescale of collapse are major con- 
tributors to the overall chemistry. 

A number of other workers in the area have reached similar 
conclusions. Garrod & Herbst (2006) specifically investigated 
methyl formate (HCOOCH3). They concluded that its forma- 
tion required the grain surface production of precursor molecules 
which were then released into the gas phase to produce methyl 
formate itself. This tallies with the similar behaviour seen here 
although Garrod & Herbst were able to produce abundances 
closer to observed values. They suggest a similar pattern for 
dimethyl ether, exactly as the Gas/Surface model shows. Later 
applications of the same model produced comparable results. 
Aikawa et al. (I2008I I consider a collapse case where again large 
organic species (or their immediate precursors) are formed on 
grain mantles before evaporation into the gas phase. Observed 
abundances can be both well modelled or over and under pro- 
duced depending on starting conditions and timescale. Garrod 
et al. ([2008 ) apply their model to more general star formation 
cases and note specifically that the longer the time available 
for grain surface chemistry to progress the more complex the 
eventual gas phase chemistry becomes. In this case agreement 
with observations is seen for abundances of ethanol and dimethyl 
ether though not for some other complex molecules. Hassel et al. 
(I2OO8I 1 specifically model L1529, a "lukewarm" corino, which 
may be the prototype of a new class of such objects. Their analy- 
sis demonstrates that, at least in the case of LI 529, there is some 
difficulty in establishing definitvely how a particular chemical 
state came to evolve and a number of scenarios are possible de- 
pending on the length of time the corino spent in each particular 
phase. 

In conclusion then grain surface processing, combined with 
the variation of physical conditions modelled here, can be 
regarded as a viable method for the formation of complex 
molecules, particularly ethanol and dimethyl ether, in the envi- 
ronment found in the vicinity of a hot corino and produce abun- 
dances comparable to those actually observed. 

Acknowledgements. Astrophysics at Queens University Belfast is supported by 
a grant from STFC. P. Hall wishes to thank Professor H. Fleischmann of Cornell 
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Table 14. Acceptable gas phase abundance ranges for other molecules 



Molecule 


Minimum 


Maximum 




Abundance 


Abundance 


H2CO 


7.0x10-^ 


8.0 xlO-" 


HCOOCH3 


1.8 xlO-* 


7.2 xlO-« 


CH3CN 


6.6x10-'" 


4.5 xlO"' 



Table 15. Other molecule abundances - gas phase 











Abundance w.r.t. H -1- H2 












H2CO 


HCOOCH3 




CH3CN 




Plateau Point 


Plateau Point 


Plateau Point 


Plateau Point 


Plateau Point 


Plateau Point 








-i-2.0xl0'yrs 




-1- 2.0 xlO' yrs 






-1-2.0x10' yrs 


Model 1 


2.95 xlO" 


-U7 


2.05 xlO""'* 


4.83 xlO"'' 


1.52 xlO""^ 


4.20 xlO" 


-m 


1.92 xlO""' 


Model 2 


4.37 xlO" 


-07 


1.77 xlQ-"* 


6.12 xlO-" 


2.19 xlO-"* 


2.44x10- 


-09 


2.85 xlO-"** 


Model 3 


3.88 xlO" 


-07 


2.23 xlO-""* 


5.57 xlO-" 


2.68 xlO-"' 


2.56 xlO- 


-09 


1.53 xlO-"' 


Model 4 


8.56 xlO" 


-07 


4.62x10-°* 


4.90 xlQ-'" 


1.33 xlO-"* 


1.64 xlO" 


-09 


1.61 xlO-"' 


Model 5 


5.32 xlO" 


-07 


4.52 xlO""^ 


3.60x10-'" 


7.18 xlO-'i 


4.05 xlO" 


-10 


5.05 xlO-'3 


Model 6 


2.00 xlO- 


-06 


1.52 xlO-"* 


1.67x10-"'' 


3.01 xlO-"' 


5.61 xlO- 


-09 


1.75 xlO-"8 


Model 7 


5.01 xlO- 


-07 


1.61 xlO-"^ 


9.06x10-'" 


5.23 xlO-" 


2.34 xlO- 


-09 


2.12 xlO-" 


Model 8 


6.70 xlO- 


-07 


2.61 xlO-"* 


2.78 xlO""* 


6.82 xl0-'2 


2.58 xlO" 


-10 


1.38 xlO"" 



Table 16. Gas/Surface model - surface phase species set 



*C 


*C2 


*C2H 


*C2H2 


*C2H3 




*C2H4 




H.C2H5O 


H.C2H5OH 


*C2H6 


♦CzHfiCO 


♦CjHO 


*C3 


*C3H 


*C3H2 


*C3H3 


*C3H4 


*C3N 


*C30 


*C4 


*C4H 


*C4H2 


*CCN 


*CCO 


*CH 


*CH2 


♦CHjCN 


H.CH2CO 


*CH2NH 


*CH20C 


♦CHjOH 


*CH3 


*CH3CHO 


H.CH3CN 


*CH30 


*CH30C 


*CH30CH 


*CH30CH2 


H.CH3OCH3 


♦CHjOH 


*CH4 


*CHOC 


*CHOH 


*CHOOH 


*CN 


*CO 


*C02 


*COC 


*COH 


*H 


*H2 


*H2C3 


*H2CN 


*H2CO 


*H20 


*H202 


N.H3C3N 


*H3C4N 


♦HCjN 


*HCN 


*HCO 


*HC00CH3 


*He 


*HNC 


*HNO 


*Mg 


*N 


*N2 


*N20 


*NH 


*NH2 


*NH, 


*NO 


*N02 


*o 


*02 


*02H 


*OCN 


*0H 



Table 17. Gas/Surface model - gas phase species set 



c 


C* 


c- 


C2 


C* 


C2H 


C2H* 


C2H2 


C2HJ 


C2H3 


C2H3 


C2H3O 


C2H3O* 


C2H4 


C2H: 


C2H4O* 


C2H5 


C2H* 


C2H5O 


C2H5O* 


C2H50H 


C2H5OH* 


C2H6 


C2H* 


C2H6CO 


C2H6CO* 


C2H6O* 


C2H6OH* 


C2H,* 


C2H7O* 


C2HO 


C2HO* 


C2HOH* 


C2NJ 


C2O* 


C3 


C* 


C3H 


C3H* 


C3H2 


C3HJ 


C3H2O+ 


C3H3 


C3H* 


C3H4 


C3H* 


C3H* 


C3H6OH+ 


C3N 


C3N* 


C30 


C3O+ 


C4 


C* 


C4H 


C4H+ 


C4H2 


C4H2 


C4H^ 


C4HJ 


C4N+ 


CCN 


CCN+ 


CCNH+ 


ceo 


CH 


CH+ 


CH2 


CHJ 


CHjCN 


CH2CN+ 


CH2CO 


CH2CO+ 


CH2NH 


CH2NHJ 


CH2OC 


CH2OH 


CH2OH+ 


CH3 


CH* 


CH3CH0 


CH3CHO+ 


CH3CN 


CH3CN+ 


CH3C0+ 


CH3O 


CH3O+ 


CH3OC 


CH3OC+ 


CH3OCH 


CH30CH+ 


CH3OCH2 


CH3OCH* 


CH3OCH3 


CH30CHJ 


CH3OH 


CH3OH* 


CH3OH* 


CH4 


CH* 


CH4N* 


CH* 


CHOC 


CHOH 


CHOH* 


CHOOH 


CHOOHJ 


CN 


CN* 


CN- 


CNC* 


CO 


CO* 


CO2 


COJ 


COC 


COM 


cooch; 


H 


H* 


H" 


H2 


HJ 


H2C3 


H2C3H* 


H2C3N+ 


H2C4N+ 


H2CN 


H2CO 


H2CO* 


H2NC+ 


H2NO+ 


H2O 


H2O* 


H2O2 




H3C3N 


H3C30+ 


H3C4N 


H3C4N* 


H3CO+ 


H,0+ 


H4C2N+ 


H4C3N+ 


H4C4N+ 


H5C2O+ 


HSC2OJ 


HC20* 


HC3N 


HC3N* 


HC3O* 


HC4N+ 


HCN 


HCN+ 


HCNH+ 


HCO 


HCO* 


HCO* 


HCOOCH3 


He 


He* 


HeH+ 


HNC 


HNCO+ 


HNO 


HNO+ 


HOC* 


Mg 


Mg* 


N 


N* 


N2 


Nt 


N2H+ 


N2O 


NCO+ 


NH 


NH+ 


NH2 


NH* 


NH3 


NHt 


NH| 


NO 


NO+ 


NO, 


NOt 





0* 


0- 


0, 


0: 


0,H 


0,H- 


OCN 


OH 


OH- 


OH- 


liLliCTR 
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Reactantl 


Reactant2 


Productl 


Product2 


Activation 
Energy £„ (K) 


1 


*C 


*CO 


*CCO 




5.0 X 10^ 


2 


*C 


*co 


*COC 




5.0 X 10^ 


3 


*C 


*H 


*CH 




0.0 


4 


*c 


*N 


*CN 




0.0 


5 


*c 


*NH 


*HNC 




0.0 


6 


*c 


*NH2 


*H 


*HNC 


0.0 


7 


*c 


*o 


*C0 




0.0 


8 


*c 


*02 


*co 


*0 


0.0 


9 


*c 


*OH 


*co 


*H 


0.0 


10 


*c 


*OH 


*COH 




0.0 


11 


H.C2H30 


*H 


*CH3CH0 




0.0 


12 


H.C2H50 


*H 


♦CjHjOH 




0.0 


13 


H.C2H0 


*H 


♦CHjCO 




0.0 


14 


*cco 


*H 


♦CjHO 




0.0 


15 


*CH 


*H 


*CH2 




0.0 


16 


*CH 


*o 


*HCO 




0.0 


17 


*CH 


*OH 


*CHOH 




0.0 


18 


*CH2 


*H 


*CH3 




0.0 


19 


*CH2 


*o 


*H2CO 




0.0 


20 


*CH2 


*0H 


*CH20H 




0.0 


21 


*CH2CO 


*H 


♦CjHjO 




0.0 


22 


H.CH2OC 


*H 


*CH30C 




0.0 


23 


*CH20H 


*H 


■nCHjOH 




0.0 


24 


*CH3 


*H 


*CH4 




0.0 


25 


*CH3 


*0 


*CH30 




0.0 


26 


*CH3 


*OH 


♦CHjOH 




0.0 


27 


H.CH3CHO 


*H 


*C2H50 




0.0 


28 


*CH30 


*H 


*CH30H 




0.0 


29 


*CH30C 


*H 


♦CHjOCH 




0.0 


30 


*CH30CH 


*H 


*CH30CH2 




0.0 


31 


♦CHjOCHj 


*H 


*CH30CH3 




0.0 


32 


*CHOC 


*H 


♦CHjOC 




0.0 


33 


*CHOH 


*H 


*CH20H 




0.0 


34 


*CN 


*H 


*HCN 




0.0 


35 


*co 


*H 


*COH 




3.0 X 10^ 


36 


*co 


*H 


*HCO 




1.0 X 103 


37 


*co 


*OH 


*C02 


*H 


3.0 X 10^ 


38 


*coc 


*H 


*CHOC 




0.0 


39 


*COH 


*H 


*CHOH 




0.0 


40 


*H 


*H 


*H2 




0.0 


41 


*H 


H.H2CO 


*H2 


*HCO 


1.85 X 10^ 


42 


*H 


*HCO 


*CO 


*H2 


0.0 


43 


*H 


*N 


*NH 




0.0 


44 


*H 


*NH 


*NH2 




0.0 


45 


*H 


*NH2 


*NH3 




0.0 


46 


*H 


*N0 


*HNO 




0.0 


47 


*H 


*o 


*OH 




0.0 


48 


*H 


*0H 


*H20 




0.0 


49 


*HCO 


*0 


*CO 


*OH 


0.0 


50 


*HCO 


*0 


*C02 


*H 


2.0 X 10' 


51 


*N 


*0 


*NO 




0.0 


52 


*N 


*OH 


*HNO 




0.0 


53 


*NH 


*o 


*HNO 




0.0 


54 


*o 


*o 


*02 




0.0 
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Table 18. Gas phase reactions used for complex molecule destiuction 





Reactantl 


Reactant2 


Product 1 


Product2 


Product3 Product4 a 




p 


y 


1 


C"^ 


C2H3O 




CO 


1.50 X 




0.00 


0.0 


2 


C"^ 


C2H5O 


C2H+ 


CO 


1.50 X 


10-09 


0.00 


0.0 


3 


C"^ 


C2HO 


C2H-" 


CO 


1.50 X 


10-09 


0.00 


0.0 


4 


C"^ 


CH2OC 


C2H+ 


HCO 


1.50 X 


10-09 


0.00 


0.0 


5 


c^ 


CH2OH 


CH+ 


CO 


1.50 X 


10-09 


0.00 


0.0 


6 


C"^ 


CH3O 


CH+ 


CO 


1.50 X 


10-09 


0.00 


0.0 


7 


c* 


CH3OC 


c 


CH+ 


CO 1.50 X 


10-09 


0.00 


0.0 


8 


C"^ 


CH3OCH 


CH 


CH+ 


CO 1.50 X 


10-09 


0.00 


0.0 


9 


c* 


CH3OCH2 


CH2 


CH+ 


CO 1.50 X 


10-09 


0.00 


0.0 


10 


c* 


CHOC 


C2H"^ 


CO 


1.50 X 


10-09 


0.00 


0.0 


11 


c* 


CHOH 


C 


CHOH^ 


7.80 X 


10-10 


0.00 


0.0 


12 


c* 


CHOH 


CH 


HCO"" 


7.80 X 


10-10 


0.00 


0.0 


13 


c* 


CHOH 


CHJ 


CO 


2.34 X 


10-09 


0.00 


0.0 


14 


c+ 


COC 


Ci 


CO 


1.50 X 


10-09 


0.00 


0.0 


15 


c+ 


COH 


ct 

2 


OH 


1.50 X 


10-09 


0.00 


0.0 


16 


C2H30 


CRPHOT 




OH 


1.30 X 


10-" 


0.00 


4307.0 


17 


C2H30 


H+ 


C2H3O+ 


H 


3.00 X 


10-09 


0.00 


0.0 


18 


C2H30 


Hi 

3 


C2H40^- 


H2 


4.85 X 


10-09 


0.00 


0.0 


19 


C2H30 


HCO* 


C2H40+ 


CO 


2.20 X 


10-09 


0.00 


0.0 


20 


C2H30 


He-" 


C2H2 


He 


OH 3.00 X 


10-09 


0.00 


0.0 


21 


C2H30 





CH3 


CO2 


1.60 X 


10-" 


0.00 


0.0 


22 


C2H30+ 


ELECTR 




OH 


1.50 X 


10-"' 


-0.50 


0.0 


23 


C2H40+ 


ELECTR 


C^H^ 


H2O 


1.50 X 


10-07 


-0.50 


0.0 


24 


C2H50 


CRPHOT 


C2H4 


OH 


1.30 X 


10-" 


0.00 


4307.0 


25 


C2H50 


H+ 


C2H5O+ 


H 


3.00 X 


10-09 


0.00 


0.0 


26 


C2H50 


at 

3 


C2H(,0+ 


H2 


4.85 X 


10-09 


0.00 


0.0 


27 


C2H50 


HCO+ 


C2H6O+ 


CO 


2.20 X 


10-09 


0.00 


0.0 


28 


C2H50 


He+ 


C2H+ 


He 


OH 3.00 X 


10-09 


0.00 


0.0 


29 


C2H50 





CH4 


CO2 


H 1.50 X 


lO-'o 


0.00 


0.0 


30 


C2H0 


CRPHOT 


C2 


OH 


1.30 X 


10-" 


0.00 


4307.0 


31 


C2HO 


H+ 


C2HO* 


H 


3.00 X 


10-09 


0.00 


0.0 


32 


C2HO 


Kt 

3 


C2HOH+ 


H2 


4.85 X 


10-09 


0.00 


0.0 


33 


C2HO 


HCO+ 


C2HOH+ 


CO 


2.20 X 


10-09 


0.00 


0.0 


34 


C2H0 


He* 


Ci 
^2 


He 


OH 3.00 X 


10-09 


0.00 


0.0 


35 


C2H0 





CO 


CO 


H 1.60 X 


10-'" 


0.00 


0.0 


36 


C2HOH+ 


ELECTR 


C2H 


OH 


1.50 X 


10-"' 


-0.50 


0.0 


37 


CH20C 


CRPHOT 


CH2 


CO 


1.30 X 


10-" 


0.00 


4307.0 


38 


CH20C 


H+ 


CH2 


HCO+ 


3.00 X 


10-09 


0.00 


0.0 


39 


CH20C 


Hi 

3 


CH4 


HCO* 


4.85 X 


10-09 


0.00 


0.0 


40 


CH20C 


HCO"^ 


CH3 


CO 


CO* 2.20 X 


10-09 


0.00 


0.0 


41 


CH20C 


He+ 


CHJ 


CO 


He 3.00 X 


10-09 


0.00 


0.0 


42 


CH20C 





CH2 


CO2 


1.60 X 


10-'" 


0.00 


0.0 


43 


CH20H 


CRPHOT 


CH2 


OH 


1.30 X 


10-" 


0.00 


4307.0 


44 


CH20H 


H+ 


CH2OH+ 


H 


3.00 X 


10-09 


0.00 


0.0 


45 


CH20H 


Ht 

3 


CH3OH+ 


H2 


4.85 X 


10-09 


0.00 


0.0 


46 


CH20H 


HCO+ 


CH3OH+ 


CO 


2.20 X 


10-09 


0.00 


0.0 


47 


CH20H 


He"- 


CHJ 


He 


OH 3.00 X 


10-09 


0.00 


0.0 


48 


CH20H 





CH3 


O2 


1.50 X 


10-" 


0.00 


0.0 


49 


CH20H 





H2CO 


OH 


1.00 X 


10-" 


0.00 


0.0 


50 


CH20H+ 


ELECTR 


CH2 


OH 


1.50 X 


10-"' 


-0.50 


0.0 


51 


CH30 


CRPHOT 


CH3 





1.30 X 


10-" 


0.00 


4307.0 


52 


CH30 


H+ 


CH3O+ 


H 


3.00 X 


10-09 


0.00 


0.0 


53 


CH30 


at 

3 


CH3OH+ 


H2 


4.85 X 


10-09 


0.00 


0.0 


54 


CH30 


HCO* 


CH3OH+ 


CO 


2.20 X 


10-09 


0.00 


0.0 


55 


CH30 


He+ 


CH+ 


He 


3.00 X 


10-09 


0.00 


0.0 


56 


CH30 





CH3 


O2 


1.50 X 


10-" 


0.00 


0.0 


57 


CH30 





H2CO 


OH 


1.00 X 


10-" 


0.00 


0.0 


58 


CH30+ 


ELECTR 


CH3 





1.50 X 


10-"' 


-0.50 


0.0 


59 


CH30C 


CRPHOT 


CH3 


CO 


1.30 X 


10-" 


0.00 


4307.0 


60 


CH30C 


H+ 


CH3OC+ 


H 


3.00 X 


10-09 


0.00 


0.0 


61 


CH30C 


H^ 


CH3OCH+ 


H2 


4.85 X 


10-09 


0.00 


0.0 


62 


CH30C 


HCO+ 


CH3OCH+ 


CO 


2.20 X 


10-09 


0.00 


0.0 


63 


CH30C 


He+ 


CH+ 


CO 


He 3.00 X 


10-09 


0.00 


0.0 


64 


CH30C 





CH3O 


CO 


1.00 X 


10-" 


0.00 


0.0 
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Table 18. Gas phase reactions used for complex molecule destiuction - continued 



15 





Reactantl 


Reactant2 


Product 1 


Product2 


Product3 


Product4 


Q: 


B 


y 


65 


CH,OC+ 


ELECTR 


CH3 


CO 






1.50 X 10"'" 


-0.50 


0.0 


66 


CH3OCH 


CRPHOT 


CH3 


HCO 






1.30 X 10"" 


0.00 


4307.0 


67 


CH3OCH 


H+ 


CH3OCH* 


H 






3.00 X 10"°' 


0.00 


0.0 


68 


CH3OCH 


HJ 

3 


CH3OCHJ 


H, 






4.85 X 10""' 


0.00 


0.0 


69 


CH3OCH 


HCO+ 


CH3OCHJ 


CO 






2.20 X 10"°' 


0.00 


0.0 


70 


CH3OCH 


He+ 


CHt 


CO 


He 




3.00 X 10"°' 


0.00 


0.0 


71 


CH3OCH 





CH3O 


CO 


H 




1.00 X 10"" 


0.00 


0.0 


72 


CH3OCH+ 


ELECTR 


CH3 


CO 


H 




1.50 X 10"°^ 


-0.50 


0.0 


73 


CH3OCH2 


CRPHOT 


CH3 


H2C0 






1.30 X 10"" 


0.00 


4307.0 


74 


CH3OCH2 


H+ 


CH3OCH* 


H 






3.00 X 10"°' 


0.00 


0.0 


75 


CH3OCH2 


HJ 


CH3OCH* 


H9 






4.85 X 10"°' 


0.00 


0.0 


76 


CH3OCH2 


HCO* 


CH3OCHJ 


CO 






2.20 X 10"°' 


0.00 


0.0 


77 


CH3OCH2 


He+ 


CH+ 


CO 


H9 


He 


3.00 X 10"°' 


0.00 


0.0 


78 


CH3OCH2 





CH3O 


H2C0 






1.00 X 10"" 


0.00 


0.0 


79 


CH3OCH2 


ELECTR 


CH3 


H2C0 






1.50 X 10"°^ 


-0.50 


0.0 


80 


CH3OCH+ 


ELECTR 


CH3 


CH30 






1.50 X 10"°^ 


-0.50 


0.0 


81 


CHOC 


CRPHOT 


CH 


CO 






1.30 X 10"" 


0.00 


4307.0 


82 


CHOC 


H+ 


CH2 


CO* 






3.00 X 10"°' 


0.00 


0.0 


83 


CHOC 


Ht 

3 


CH3 


HCO* 






4.85 X 10"°' 


0.00 


0.0 


84 


CHOC 


HCO* 


CH2 


CO 


CO* 




2.20 X 10"°' 


0.00 


0.0 


85 


CHOC 


He* 


CH* 


CO 


He 




3.00 X 10"°' 


0.00 


0.0 


86 


CHOC 





CH 


C02 






1.60 X 10"'° 


0.00 


0.0 


87 


CHOH 


CRPHOT 


CO 


H2 






1.30 X 10"" 


0.00 


2659.0 


88 


CHOH 


H+ 


CHOH* 


H 






2.96 X 10"°' 


0.00 


0.0 


89 


CHOH 


H+ 


CO* 


H 


H, 




1.06 X 10"°' 


0.00 


0.0 


90 


CHOH 


H+ 


Ho 


HCO* 






3.57 X 10"°' 


0.00 


0.0 


91 


CHOH 


HJ 


H9 


H3CO* 






6.30 X 10"°' 


0.00 


0.0 


92 


CHOH 


HCO+ 


CO 


H3CO+ 






3.30 X 10"°' 


0.00 


0.0 


93 


CHOH 


He+ 


CH* 


He 







1.71 X 10"°' 


0.00 


0.0 


94 


CHOH 


He+ 


CHOH* 


He 






9.69 X 10"'° 


0.00 


0.0 


95 


CHOH 


He* 


CO* 


H2 


He 




1.88 X 10"°' 


0.00 


0.0 


96 


CHOH 


He* 


H 


HCO* 


He 




1.14 X 10"°' 


0.00 


0.0 


97 


CHOH 





HCO 


OH 






1.78 X 10"" 


0.57 


1390.0 


98 


CHOH+ 


ELECTR 


CO 


H 


H 




5.00 X 10"°'' 


-0.50 


0.0 


99 


CHOH+ 


ELECTR 


H 


HCO 






1.00 X 10"°'' 


-0.50 


0.0 


100 


CHOH+ 


ELECTR 


H2CO 


PHOTON 






1.10 X 10"'° 


-0.70 


0.0 


101 


COC 


CRPHOT 


c 


CO 






1.30 X 10"'^ 


0.00 


4307.0 


102 


COC 


H+ 


CH 


CO* 






3.00 X 10"°' 


0.00 


0.0 


103 


COC 


"3 


CH3 


CO* 






4.85 X 10"°' 


0.00 


0.0 


104 


COC 


HCO* 


CH 


CO 


CO* 




2.20 X 10"°' 


0.00 


0.0 


105 


COC 


He* 


C* 


CO 


He 




3.00 X 10"°' 


0.00 


0.0 


106 


COC 





CO 


CO 






1.60 X 10"'° 


0.00 


0.0 


107 


COH 


CRPHOT 


C 


OH 






1.30 X 10"'^ 


0.00 


4307.0 


108 


COH 


H+ 


CH 


OH* 






3.00 X 10"°' 


0.00 


0.0 


109 


COH 


H* 


CH3 


OH* 






4.85 X 10"°' 


0.00 


0.0 


110 


COH 


HCO* 


CH 


CO* 


OH 




2.20 X 10"°' 


0.00 


0.0 


111 


COH 


He* 


C* 


He 


OH 




3.00 X 10"°' 


0.00 


0.0 


112 


COH 





CO 


OH 






1.60 X 10"'° 


0.00 


0.0 
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Table 19. Model 1 - Free fall collapse from 100 cm"' - most abundant gas phase species 



Plateau Point Plateau Point + 2.0 xlO^ years 
4.389 X 10^ years 4.589 x 10^ years 





Species 


41 1 / *TT TT\ 

Abundance (w.r.t H + H2) 




Species 


Abundance (w.r.t H + H2) 


1 


Hz 


5.000 X 10""' 


1 


H2 


5.000 X 10""' 


2 


He 


1.401 X 10""' 


2 


He 


1.401 X 10""' 


3 


H2O 


1.554 X lO"** 


3 


H2O 


1.447 X 10""^ 


4 


CH4 


5.260 X 10""^ 


4 


CH4 


3.980 X 10""^ 


5 


CH3OH 


1.116 X 10""^ 


5 


CO2 


1.462 X 10""^ 


6 


NH3 


8.369 X 10""* 


6 


N2 


7.593 X 10""* 


7 


N2 


6.114 X 10""* 


7 


HCN 


2.847 X 10""* 


8 


CO 


4.929 X 10""* 


8 


HNC 


1.607 X 10""* 


9 


CO2 


1.352 X 10""* 


9 


C4H 


8.936 X 10""'' 


10 


NH2 


3.999 X 10""'' 


10 


CH2CN 


8.380 X 10""'' 


11 


CH3 


3.254 X lO""' 


11 


CH3OCH3 


8.361 X 10""'' 


12 


CH3OCH3 


3.006 X 10""'' 


12 


C3H3 


6.461 X lO""' 


13 


H2CO 


2.946 X 10""^ 


13 


HC3N 


5.801 X 10""'' 


14 


O2 


2.416 X 10""'' 


14 


CH3OH 


4.315 X 10""'' 


15 


HCN 


1.916 X 10""'' 


15 


CiH 


3.657 X 10""^ 


16 


C4H 


1.545 X 10""^ 


16 


C2 


2.372 X 10""^ 


17 


HNC 


1.193 X 10""^ 


17 


CH3CN 


1.923 X 10""^ 


18 


C3H 


8.901 X 10""* 


18 


C4H2 


6.160 X 10""* 


19 


CH 


4.796 X 10""» 


19 


CH 


5.320 X 10""* 


20 


NH 


2.391 X 10""* 


20 


C3H2 


5.081 X 10""* 


21 


C3H3 


1.873 X 10""* 


21 


C3N 


3.511 X 10""* 


22 


H 


1.364 X 10""* 


22 


C3 


3.456 X 10""* 


23 


CH2 


9.769 X 10""" 


23 


CH3 


2.655 X 10""* 


24 


C2H4 


8.988 X 10""" 


24 


CCN 


2.442 X 10""* 



Table 20. Model 2 - Free fall collapse from 1,000 cm ^ - most abundant gas phase species 

Plateau Point Plateau Point + 2.0 xlO'' years 
1.388 X 10'' years 1.588 x 10'' years 



Species Abundance (w.r.t H + H2) Species Abundance (w.r.t H + H2) 



1 


H2 


5.000 X 10""' 


1 


H2 


5.000 X 10""' 


2 


He 


1.401 X 10""' 


2 


He 


1.401 X 10""' 


3 


H2O 


1.375 X lO""* 


3 


H2O 


1.300 X 10""^ 


4 


CH4 


3.880 X 10""^ 


4 


CH4 


3.118 X 10""^ 


5 


CH3OH 


1.834 X 10""5 


5 


CO2 


2.121 X 10""^ 


6 


CO 


8.797 X 10""* 


6 


N2 


8.945 X 10""* 


7 


N2 


8.063 X 10""* 


7 


CH3 


3.178 X 10""* 


8 


NH3 


4.633 X 10""* 


8 


H2CO 


1.768 X 10""* 


9 


CO2 


3.475 X 10""* 


9 


C4H 


1.450 X 10""* 


10 


O2 


1.413 X 10""* 


10 


HCN 


1.444 X 10""* 


11 


CH3 


4.982 X 10""^ 


11 


HNC 


1.382 X 10""* 


12 


H2CO 


4.366 X 10""' 


12 


CH3OCH3 


5.813 X 10""' 


13 


CH3OCH3 


4.110 X 10""' 


13 


C2 


4.082 X 10""' 


14 


NH2 


3.131 X 10""' 


14 


CH3OH 


3.799 X 10""' 


15 


C4H 


2.115 X lO""' 


15 


C3H 


3.768 X 10""' 


16 


HCN 


1.361 X 10""' 


16 


CH,CN 


3.276 X 10""' 


17 


HNC 


8.830 X 10""* 


17 


C3H3 


2.841 X 10""' 


18 


C3H 


7.771 X 10""* 


18 


HC3N 


2.025 X 10""' 


19 


CH 


3.446 X 10""* 


19 


C3H2 


8.748 X 10""* 


20 


NH 


2.147 X 10""* 


20 


C3 


8.348 X 10""* 


21 


C2H4 


1.672 X 10""* 


21 


CH 


6.791 X 10""* 


22 


H 


1.329 X 10""* 


22 


C4H2 


6.266 X 10""* 


23 


C2H2 


1.205 X 10""* 


23 


CO 


4.991 X 10""* 


24 


OCN 


1.130 X 10""* 


24 


C2H50H 


3.318 X 10""* 
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Table 21. Model 3 - Free fall collapse from 10,000 cm"' - most abundant gas phase species 



Plateau Point Plateau Point + 2.0 xlO^ years 
4.381 X 10^ years 6.381 x 10^ years 





Species 


jil 1 / *TT TT\ 

Abundance (w.r.t H + H2) 




Species 


Abundance (w.r.t H + H2) 


1 


Hz 


5.000 X 10""' 


1 


H2 


5.000 X 10""' 


2 


He 


1.401 X 10""' 


2 


He 


1.401 X 10""' 


3 


H2O 


1.494 X lO"** 


3 


H2O 


1.403 X 10""^ 


4 


CH4 


4.801 X 10""' 


4 


CH4 


3.849 X 10""' 


5 


CH3OH 


1.555 X 10""' 


5 


CO2 


1.667 X 10""' 


6 


N2 


6.961 X 10""* 


6 


N2 


8.063 X 10""* 


7 


NH3 


6.807 X lO""*" 


7 


HCN 


2.385 X 10""* 


8 


CO 


4.811 X 10"""^ 


8 


HNC 


1.355 X 10""* 


9 


CO2 


1.700 X 10""'^ 


9 


CH3OCH3 


1.038 X 10""* 


10 


O2 


6.601 X 10""'' 


10 


C4H 


8.916 X 10""' 


11 


CH3 


4.107 X 10""' 


11 


CH2CN 


6.885 X 10""' 


12 


H2CO 


3.875 X 10""'' 


12 


C3H3 


5.968 X lO""' 


13 


CH3OCH3 


3.630 X 10""^ 


13 


HC3N 


5.436 X 10""' 


14 


NH2 


3.204 X 10""'' 


14 


CH3OH 


5.232 X 10""' 


15 


HCN 


1.606 X 10""'' 


15 


C3H 


3.802 X 10""' 


16 


C4H 


1.090 X 10""' 


16 


C2 


2.579 X 10""' 


17 


HNC 


8.841 X 10""** 


17 


CH3CN 


1.528 X 10""' 


18 


C3H 


6.344 X 10""* 


18 


C4H2 


6.469 X 10""* 


19 


CH 


4.268 X 10""* 


19 


C3H2 


5.613 X 10""* 


20 


NH 


2.196 X 10""* 


20 


CH 


5.452 X 10""* 


21 


OCN 


1.872 X 10""* 


21 


C3 


3.771 X 10""* 


22 


C3H3 


1.502 X 10""* 


22 


CH3 


3.374 X 10""* 


23 


H 


1.335 X 10""* 


23 


C3N 


3.340 X 10""* 


24 


C2H4 


1.332 X 10""* 


24 


H2CO 


2.230 X 10""* 



Table 22. Model 4 - Retarded collapse from 100 cm ' - most abundant gas phase species 



Plateau Point Plateau Point + 2.0 xlO'' years 



2.193 X 10' years 




2.21 


13 X l(y yeans 






Species 


Abundance (w.r.t H + H2) 




Species 


Abundance (w.r.t H + H2) 


1 


H2 


5.000 X 10""' 


1 


H2 


5.000 X 10""' 


2 


He 


1.401 X 10""' 


2 


He 


1.400 X 10""' 


3 


H2O 


1.158 X 10"** 


3 


H2O 


1.021 X 10""^ 


4 


CH4 


3.190 X 10""' 


4 


CO2 


2.822 X 10""' 


5 


CO 


1.416 X 10""' 


5 


CH4 


2.185 X 10""' 


6 


CO2 


1.155 X 10""' 


6 


CO 


1.184 X 10""' 


7 


N2 


1.033 X 10""' 


7 


N2 


1.034 X 10""' 


8 


CH3OH 


8.562 X 10""* 


8 


H2CO 


4.623 X 10""* 


9 


O2 


5.938 X 10""* 


9 


C4H 


9.170 X 10""' 


10 


CH3OCH3 


8.609 X 10""' 


10 


HCN 


2.600 X 10""' 


11 


H2CO 


8.562 X 10""' 


11 


HNC 


2.474 X 10""' 


12 


CH3 


6.489 X 10""' 


12 


C3H 


1.636 X 10""' 


13 


C4H 


5.230 X 10""' 


13 


C3H3 


1.141 X 10""' 


14 


NH3 


2.774 X 10""' 


14 


CH3OH 


7.979 X 10""* 


15 


C3H 


1.241 X 10""' 


15 


C3H2 


7.747 X 10""* 


16 


HCN 


1.232 X 10""' 


16 


CH3OCH3 


6.864 X 10""* 


17 


HNC 


1.008 X 10""' 


17 


CH2CN 


6.811 X 10""* 


18 


C2H2 


8.272 X 10""* 


18 


C2H5OH 


6.214 X 10""* 


19 


NH2 


7.967 X 10""* 


19 


HC3N 


4.022 X 10""* 


20 


OCN 


5.908 X 10""* 


20 


C4H2 


2.109 X 10""* 


21 


C3H2 


4.136 X 10""* 


21 


C2 


1.836 X 10""* 


22 


C2H4 


3.294 X 10""* 


22 


CH 


1.541 X 10""* 


23 


CH 


2.817 X 10""* 


23 


C3 


1.504 X 10""* 


24 


C2 


2.198 X 10""* 


24 


H 


1.501 X 10""* 
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Table 23. Model 5 - Retarded collapse from 1,000 cm"' - most abundant gas phase species 



Plateau Point Plateau Point + 2.0 xlO^ years 
6.930 X 10^ years 7.130 x 10^ years 





Species 


jil 1 / *TT TT\ 

Abundance (w.r.t H + H2) 




Species 


Abundance (w.r.t H + H2) 


1 


Hz 


5.000 X 10""' 


1 


H2 


5.000 X 10""' 


2 


He 


1.400 X 10""' 


2 


He 


1.400 X 10""' 


3 


HjO 


6.859 X 10""' 


3 


H2O 


4.986 X 10""' 


4 


CO 


4.396 X 10""' 


4 


CO2 


4.113 X 10""' 


5 


CO2 


1.946 X 10""' 


5 


CO 


2.959 X 10""' 


6 


N2 


1.056 X 10""' 


6 


N2 


1.064 X 10""' 


7 


O2 


8.515 X lO""*" 


7 


O2 


6.551 X 10""*" 


8 


CH3OH 


5.436 X lO""*^ 


8 


CH4 


1.308 X 10""*^ 


9 


CH3OCH3 


7.695 X 10""'' 


9 


H2CO 


4.522 X 10""^ 


10 


CH3 


5.394 X 10""'' 


10 


C2H4 


6.840 X 10""* 


11 


H2CO 


5.316 X lO""' 


11 


C2H2 


2.861 X 10""* 


12 


CH4 


4.808 X 10""'' 


12 





2.465 X 10""* 


13 


C2H2 


9.598 X 10""* 


13 


CHOOH 


1.970 X 10""* 


14 


C4H 


6.895 X 10""* 


14 


H 


1.893 X 10""* 


15 


HCN 


4.602 X 10""* 


15 


CH2CN 


9.355 X 10""" 


16 


NH3 


4.573 X 10""* 


16 


Mg 


6.969 X 10""" 


17 


OCN 


3.380 X 10""* 


17 


C4H 


5.114 X 10""" 


18 


HNC 


3.317 X 10""* 


18 


OH 


3.563 X 10""" 


19 


C3H 


2.731 X 10""* 


19 


HCN 


3.234 X 10""" 


20 


NH2 


2.527 X 10""* 


20 


HNC 


2.975 X 10""" 


21 


C3H2 


2.332 X 10""* 


21 


N 


2.669 X 10""" 


22 


H 


1.423 X 10""* 


22 


NH 


2.406 X 10""" 


23 


Mg 


6.970 X 10""" 


23 


NH3 


1.498 X 10""" 


24 


NH 


5.941 X 10""" 


24 


CH3OH 


1.343 X 10""" 



Table 24. Model 6 - Retarded collapse from 10,000 cm ^ - most abundant gas phase species 



Plateau Point 

2.19 8 X 10^ years 

Species Abundance (w.r.t H + H2) 



Plateau Point + 2.0 xlO^ years 
2.398 X 10'' years ' 



Species 



Abundance (w.r.t H + H2) 



1 


H2 


5.000 X 10""' 


1 


H2 


5.000 X 10""' 


2 


He 


1.401 X 10""' 


2 


He 


1.401 X 10""' 


3 


H2O 


1.391 X 10"** 


3 


H2O 


1.258 X lO"*' 


4 


CH4 


4.169 X 10""' 


4 


CH4 


3.347 X 10""' 


5 


CH3OH 


1.094 X 10""' 


5 


CO2 


2.455 X 10""' 


6 


N2 


1.007 X 10""' 


6 


N2 


9.940 X 10""* 


7 


CO2 


7.513 X 10""'' 


7 


C4H 


2.133 X 10""* 


8 


CO 


2.971 X 10""* 


8 


HCN 


6.115 X 10""^ 


9 


H2CO 


1.995 X 10""'^ 


9 


C3H 


4.118 X 10""'' 


10 


O2 


1.747 X 10""'^ 


10 


HNC 


3.100 X 10""' 


11 


CH3OCH3 


1.639 X 10""'^ 


11 


HC3N 


2.558 X 10""' 


12 


CH3 


1.302 X 10""* 


12 


C3H3 


2.388 X 10""' 


13 


NH3 


5.736 X 10""^ 


13 


C2 


2.148 X 10""' 


14 


C4H 


3.498 X 10""^ 


14 


CH2CN 


2.026 X 10""' 


15 


HCN 


2.723 X 10""^ 


15 


CH3OCH3 


1.920 X 10""' 


16 


HNC 


1.553 X 10""^ 


16 


CH3OH 


1.656 X 10""' 


17 


C3H 


1.509 X 10""^ 


17 


C3H2 


9.933 X 10""* 


18 


OCN 


7.426 X 10""* 


18 


C4H2 


9.856 X 10""* 


19 


C3H3 


7.264 X 10""* 


19 


C3 


5.645 X 10""* 


20 


NH2 


6.593 X 10""* 


20 


CH 


5.217 X 10""* 


21 


C3H2 


4.091 X 10""* 


21 


C3N 


2.139 X 10""* 


22 


CH 


4.007 X 10""* 


22 


CH3CN 


1.751 X 10""* 


23 


C2H4 


3.178 X 10""* 


23 


CH3 


1.596 X 10""* 


24 


C2 


2.231 X 10""* 


24 


H2CO 


1.518 X 10""* 
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Table 25. Model 7 - Free fall collapse from 10,000 cm"' - H2 starting abundance - most abundant gas phase species 



Plateau Point Plateau Point + 2.0 xlO^ years 
4.381 X 10^ years 6.381 x 10^ years 





Species 


jil 1 / *TT TT\ 

Abundance (w.r.t H + H2) 




Species 


Abundance (w.r.t H + H2) 


1 


H2 


5.000 X 10""' 


1 


H2 


5.000 X 10""' 


2 


He 


1.400 X 10""' 


2 


He 


1.400 X 10""' 


3 


H2O 


8.079 X 10""' 


3 


H2O 


5.528 X 10""' 


4 


CO 


3.734 X 10""' 


4 


CO2 


4.766 X 10""' 


5 


CO2 


2.562 X 10""' 


5 


CO 


2.288 X 10""' 


6 


N2 


1.054 X 10""' 


6 


N2 


1.064 X 10""' 


7 


CH4 


2.491 X lO""*" 


7 


CH4 


1.801 X 10""*" 


8 


O2 


2.137 X lO""*^ 


8 


O2 


8.088 X 10""'' 


9 


CH3OH 


8.713 X 10""'' 


9 


H2CO 


1.609 X 10""'' 


10 


C4H 


7.661 X 10""'' 


10 





2.585 X 10""* 


11 


C3H 


5.091 X lO""' 


11 


C2H2 


2.445 X 10""* 


12 


H2CO 


5.005 X 10""'' 


12 


C4H 


2.276 X 10""* 


13 


CH3OCH3 


1.936 X lO""' 


13 


CHOOH 


2.094 X 10""* 


14 


C2H2 


1.038 X 10""'' 


14 


CH2CN 


2.088 X 10""* 


15 


CCN 


5.875 X 10""^ 


15 


H 


1.890 X 10""* 


16 


CH3 


5.190 X lO""** 


16 


C2H4 


7.726 X 10""' 


17 


HCN 


4.327 X lO""** 


17 


Mg 


6.969 X 10""' 


18 


C3 


3.276 X 10-"* 


18 


C3H 


4.168 X 10""' 


19 


NH3 


2.411 X 10""* 


19 


OH 


3.816 X 10-"' 


20 


C2 


2.391 X 10""* 


20 


CH3OH 


2.722 X 10-"' 


21 


C3H2 


2.371 X 10""* 


21 


N 


2.559 X 10""' 


22 


HNC 


2.314 X 10""* 


22 


NH 


2.388 X 10""' 


23 


C2H5OH 


2.301 X 10""* 


23 


NH3 


1.793 X 10""' 


24 


HC3N 


1.915 X 10""* 


24 


HCN 


1.226 X 10""' 



Table 26. Model 8 - Retarded collapse from 10,000 cm ^ - H2 starting abundance - most abundant gas phase species 

Plateau Point Plateau Point + 2.0 xlO'' years 
2.198 X 10'' years 2.398 x 10'' yeans 



Species Abundance (w.r.t H + H2) Species Abundance (w.r.t H + H2) 



1 


H2 


5.000 X 10""' 


1 


H2 


5.000 X 10""' 


2 


He 


1.400 X 10""' 


2 


He 


1.400 X 10""' 


3 


H2O 


6.243 X 10""' 


3 


CO2 


5.294 X 10""' 


4 


CO 


3.772 X 10""' 


4 


H26 


3.453 X 10""' 


5 


CO2 


2.777 X 10-"' 


5 


CO 


1.839 X 10-"' 


6 


N2 


1.057 X 10""' 


6 


N2 


1.065 X 10""' 


7 


O2 


9.298 X 10""« 


7 


O2 


8.183 X 10-"« 


8 


CH4 


1.577 X 10""« 


8 


CH4 


1.207 X 10-"« 


9 


C4H 


6.777 X 10""^ 


9 


CHOOH 


5.060 X 10""* 


10 


H2CO 


6.702 X 10""^ 


10 





3.969 X 10""* 


11 


C2H2 


2.364 X 10""^ 


11 


C2H4 


3.818 X 10""* 


12 


C3H 


2.304 X 10-"' 


12 


H2CO 


2.609 X 10-"* 


13 


CH3OH 


1.507 X 10""' 


13 


H 


1.888 X 10""* 


14 


CH3OCH3 


1.266 X 10""' 


14 


C2H2 


1.256 X 10""* 


15 


C3H2 


8.157 X 10""* 


15 


CH2CN 


1.058 X 10""* 


16 


HCN 


3.963 X 10""* 


16 


Mg 


6.968 X 10""' 


17 


CCN 


3.628 X 10""* 


17 


OH 


3.645 X 10""' 


18 


HNC 


3.490 X 10""* 


18 


N 


3.416 X 10""' 


19 


C2 


1.835 X 10-"* 


19 


CH3OH 


1.968 X 10""' 


20 


OCN 


1.683 X 10""* 


20 


NH 


1.784 X 10""' 


21 


C4H2 


1.657 X 10""* 


21 


NH3 


1.611 X 10-"' 


22 


HC3N 


1.641 X 10""* 


22 


NO 


1.345 X 10""' 


23 


H 


1.619 X 10""* 


23 


C4H 


5.672 X 10"'" 


24 


C3 


1.461 X 10""* 


24 


HNO 


4.566 X 10"'" 
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Fig. 3. Model 3 - Free fall collapse from 10,000 cm ^ - gas phase 
data 
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Fig. 5. Model 5 - Retarded collapse from 1,000 cm ^ - gas phase 
data 
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Fig. 6. Model 6 - Retarded collapse from 10,000 cm ^ 
phase data 
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Fig. 7. Model 7 - Free fall collapse from 10,000 cm-^ - H2 start- 
ing abundance - gas phase data 
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ing abundance - gas phase data 
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